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Abstract 

This report describes a system which maintains canonical expressions for 
designators under a set of equalities. Substitution is used to maintain all 
knowledge in terms of these canonical expressions. A partial order on 
designators, termed the better-name relation, is used in the choice of canonical 
expressions. It is shown that with an appropriate better-name relation an 
important engineering reasoning technique, propagation of constraints, can be 
implemented as a special case of this substitution process. Special purpose 
algebraic simplification procedures are embedded such that they interact 
effectively with the equality system. An electrical circuit analysis system is 
developed which relies upon constraint propagation and algebraic simplification as 
primary reasoning techniques. The reasoning is guided by a better-name relation 
in which referentially transparent terms are preferred to referentially opaque 
ones. Multiple description of subcircuits are shown to interact strongly with the 
reasoning mechanisms. 

A conceptual distinction is made between propositional deduction and 
the instantiation of quantified knowledge. A special purpose system is used for 
the simplified domain of propositional deduction which incorporates many useful 
features. This system, termed a truth maintenance system or TMS, keeps track 
of justifications for deductions which are used to generate explanations of 
deduced beliefs. The TMS also handles the retraction of assumptions when 
contradictions arise. When assumptions are retracted all deductions which 
depended on them are retracted in an efficient incremental manner. Assumptions 
that certain algebraic quantities are non-zero play an important role in reasoning 
about algebraic constraints. 
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Overview 

In almost all reasoning systems, some method of handling the equality of 
objects must be devised. The problems associated with equality, or multiple 
expressions which have the same referent, are usually treated as an orthogonal 
issue to other problem solving or reasoning techniques. The approach taken here 
is that multiple descriptions for objects and the issues surrounding equality are of 
central importance in problem solving itself. A uniform algorithm for handling 
equality is developed here which assigns canonical names to equivalence classes of 
expressions. Specific reasoning techniques, such as propagation of constraints, are 
built up using this algorithm as a primary deductive mechanism. Propagation of 
constraints is then used as the major reasoning mechanism in an electronic 
analysis system. Figure 1 shows an overall breakdown of the reasoning system. 
The remainder of this section describes the various systems and the way they 
depend on one another. 
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Figure 1. The Major Subsystems. 

Underlying all the systems described here is a truth maintenance system 
(TMS) which handles all propositional deduction [Doyle 79] [McAUester 78]. The 
TMS keeps track of justifications for all deduced truth values. These 
justifications are used to generate explanations for any deduced truth value for a 
proposition. Such explanations give the user a better understanding of the 
deductions being performed and therefore more confidence in the results. The 
justifications are also used to incrementally update the truth value of all 
propositions when assumptions are added or retracted. When contradictions 
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occur the assumptions underlying them can be pinpointed and the negation of 
one of these assumptions can be deduced in a process termed dependency 
directed backtracking [Stallman & Sussman 77]. All of these functions, 
explanation veneration, incremental modification, and dependency directed 
backtracking are performed by special purpose algorithms in the TMS. 

All of the TMS functions described above were incorporated into the 
original version of a TMS formulated by Jon Doyle [Doyle 78]. In addition to 
these functions the TMS is here viewed as an important deductive component of 
the overall system. A distinction is made between propositional deduction and 
the instantiation of quantified knowledge. Quantified knowledge is taken to be 
any knowledge which can be usefully applied to a large number of specific 
individuals (knowledge concerning the existence of certain types of individuals is 
not dealt with). The instantiation of quantified knowledge is the act of 
generating specific knowledge about an individual corresponding to that 
quantified knowledge. Often instantiation can be made simply in the hope that 
the resulting propositional knowledge will be useful in propositional deduction. 
Since the propositional reasoning mechanisms in the TMS are rather 
straightforward the major control issue is the control of the instantiation of 
quantified knowledge. 

A restriction on the current version of the systems is that it is not 
possible to add quantified knowledge dynamically during reasoning and have it 
retrospectively instantiated with all the appropriate items. None of the 
applications described here require such an ability, but there is no theoretical 
reason that this limitation could not be overcome within a very similar 
framework to the one discussed here. 

The equality system described in chapter one deals with propositions 
asserting equalities between "designators". My use of the word "designator" is 
equivalent to the use of the word "term" in standard systems of formal logic. 
The word "designator" is used here for no better reason than that it seems to 
focus ones attention on the distinction between a designator and its referent, and 
therefore seems more natural in a system centered around equality and the use of 
multiple designators. Each designator is maintained as an independent entity 
even though an equality may state that it has the same referent as some other 
designator. The equalities in the system determine equivalence classes of 
designators. These classes are incrementally maintained as equalities are added 
and removed from the system. The equality system maintains a canonical name 
for each equivalence class which is one of the designators in that class. At any 
point the user of the system can ask for the canonical name of some designators 
equivalence class via a "what-is" function. A substitution process is also 
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controlled by the equality system which essentially insures that all knowledge is 
represented in terms of canonical names. 

The equality system can be viewed as an algorithm for controlling the 
instantiation of a specific body of quantified knowledge. The quantified 
knowledge involved is the transitivity of equality and the validity of the 
substitution of equals for equals. The equality system instantiates this knowledge 
by creating new equalities and telling the TMS that the new equalities are 
implied by other equalities already in the system. Of course the mechanism for 
controlling this instantiation process is largely determined by the set of equalities 
believed by the system. An important observation about this particular control 
of instantiation is that it is object oriented. Thus, while the system is based on 
an assertional data base, it simulates the behavior of systems built on object 
oriented data structures. 

A major reasoning technique employed by electrical engineers is 
propagation of constraints [Steele & Sussman 78]. It is shown in chapter one that 
this process can be implemented as a special case of the algorithms for handling 
equality. The basic step in constraint propagation is a shift in the canonical 
name of a class induced by a substitution of terms into some designator in that 
class. The shift in canonical name can then induce further substitutions thus 
propagating values around a constraint network. 

There are several mechanisms which have been developed to allow the 
procedural embedding of knowledge. The most straightforward example of this is 
the use of arithmetic operators. For example consider the designator (+ 1 2). 
The arithmetic operator + can be applied to 1 and 2 to give a new numerical 
designator, 3, which has the same referent as the sum. Thus equalities can be 
generated by instantiating knowledge about the functions used in designators and 
this knowledge can be contained in procedures attached to these functions. 
Algebraic simplification is done essentially the same way. A few other 
mechanisms for the procedural embedding and special purpose controlled 
instantiation within the equality system are discussed in chapter two. 

Chapter three introduces electronic analysis. Knowledge about simple 
electronic circuits is naturally embedded in TMS predicates and propagation of 
constraints provides the primary reasoning mechanism. Canonical names are 
chosen using a "better-name" relation on designators and several important points 
are made about the structure of the better-name relation necessary to produce a 

transparent designators must be better names than referentially opaque ones. 
Another important observation is that arbitrary quantities (ones that stand for an 
unspecified and arbitrary value such as a ground potential or an input voltage) 
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should be better names than referentially opaque designators which are not 
arbitrary (arbitrary values are not used in this system to derive quantified 
knowledge). 

A mechanism for using multiple descriptions of parts of electronic 
circuits, termed "slices", has been proposed by Gerald Sussman [Sussman 77]. 
This mechanism has been implemented in the equality system via the creation of 
equivalences between the terminals of equivalent circuits. The introduction of 
such equivalences fits naturally into the equality framework. A restructuring of 
constraints which corresponds to redrawing circuit topologies is also shown to be 
important for propagation of constraints and the effectiveness of slices. These 
techniques are discussed in chapter three. 

Chapter four gives some of the algorithmic details which were glossed 
over in the previous chapters. A description of some related previous work is 
also presented. 

Everything described here has been fully implemented and is currently 
(July 79) running on both the PDP-10 and the LISP machines at the MIT AI 
Lab. 
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Chapter I 
The TMS, Canonical Naming, and Constraints 



Using the TMS 

The truth maintenance system or TMS, is a special purpose system for 
handling propositional deduction [Doyle 79] [McAllester 78]. The TMS used here 
should be thought of as a system for enforcing a set of logical relations on a set 
of propositions. Each proposition is represented in the TMS by a TMS node. 
Each such node can be in one of three possible truth states, true, false, or 
unknown. The logical relations in the TMS should be thought of as constraints 
on the truth values associated with the nodes in that relation. Each relation 
(constraint) is represented internally as a disjunction of terms such as (or (not P) 
(not 0) R). The details of the TMS algorithms are discussed further in chapter 
four. 

When a deduction is made within the TMS a justification for that 
deduction is recorded. The resulting justifications are used to generate 
explanations for any deduced assertion when such explanations are requested by 
the user. These explanations give the user insight into the way a deduction was 
made and therefore greater confidence in the results. The justifications are far 
more important than this however. They allow the TMS to incrementally modify 
the set of deduced truth values when assumptions are retracted. When 
contradictions occur in the system the assumptions underlying that contradiction 
can be pinpointed via the justifications, and the negations of one these 
assumptions can be deduced to remove the contradiction. This controlled 
backing out from a contradiction has been termed dependency directed 
backtracking [Stallman & Sussman 77]. 

With the introduction of a special purpose propositional deduction 
system (the TMS) a sharp distinction is being made between propositional 
deduction and the instantiation of quantified knowledge. Quantified knowledge is 
taken to be any knowledge which can be usefully applied to a large number of 
specific individuals (knowledge concerning the existence of certain types of 
individuals is not dealt with). The instantiation of quantified knowledge is the 
act of generating specific knowledge about an individual corresponding to that 
quantified knowledge. In the present system all such specific knowledge is 
represented as either a logical relation among TMS nodes or as the assignment of 
a truth value to such a node. For example consider the quantified knowledge 
represented by the sentence: Vx(mammal(x) -> warm-blooded(x)). This could be 
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instantiated with "Fred" to give the propositional relation: mammal(Fred) -> 
warm-blooded(Fred). Now if the TMS believes that the proposition 
mammal(Fred) is true (the node corresponding to this assertion has a truth value 
of "true"), it will deduce that warm-blooded(Fred) is also true. However if the 
TMS does not believe that Fred is a mammal, i.e. either it believes that Fred is 
not a mammal or it has not committed itself one way or the other, then it will 
simply remember the implication and use it whenever it can. The control of 
instantiation is the major control issue in the reasoning system. 

The simplest embodiment of quantified knowledge is in "TMS 
predicates". These predicates are lisp functions which return a TMS node 
representing a proposition. For example a mammal predicate can be defined as a 
function of one variable which returns the node representing the proposition that 
the argument passed is a mammal. Quantified knowledge is embodied in such 
predicates in that instantiations are done as side effects the first time a predicate 
is applied to some set of arguments. Thus the example of quantified knowledge 
above which states that all mammals are warm blooded could be attached to the 
mammal predicate and instantiated with every object to which the mammal 
predicate is applied. Notice that this is not "antecedent deduction" in the 
classical sense as no actual deduction need be involved at all. I will call this type 
of instantiation "reference instantiation" since the quantified knowledge is 
instantiated at the point reference is made to a certain entity (such as a 
proposition representing the application of the mammal predicate). A special 
form has been defined for the creation of TMS predicates and some examples of 
TMS predicate definitions are given below: 

(defpred vertebrate (x)) 

(defpred flys (x)) 

(defpred bird (x) 

(vertebrate x) 
(flys x)) 

The body of a defpred must be a list of forms which return TMS nodes 
when evaluated. The TMS nodes thus created will be implied by the truth of the 
node returned by an application of the predicate being defined. Thus (bird 
'John) returns a TMS node which implies (vertebrate 'John) and (flys 'John). 

There are various TMS predicates which are supplied as primitives. The 
predor and predand predicates take any number of TMS nodes and return a node 
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which has been appropriately constrained in the TMS. The -> predicate takes 
two nodes and returns a TMS node which implies an implication relation between 
them. Prednot returns a node which is constrained to be equivalent to the 
negation of the node it is passed. Finally presumably returns a node representing 
the assertion that the node it was passed should be assumed to be true. This 
interpretation is implemented by using the -> predicate to create an assertion of 
the form: (-> (presumably P) P). The node representing this implication is made 
true as an assumption which can be retracted if contradictions arise. To see how 
these primitive predicates can be used consider the following definition of a 
mammal predicate: 

(defpred hairy (x)) 

(defpred female (x)) 

(defpred bears-l i ve-young (x)) 

(defpred mammal (x) 

(vertebrate x) 
(presumably (hairy x)) 
(presumably (-> (female x) 

(bears-l i ve-young x)))) 

There is a TMS primitive for setting the truth value of TMS nodes 
which can be used to declare that some object is a mammal as follows: 



or equivalently 



(set-truth (mammal 'Joe) 'true 'premise) 



(assert (mammal 'joe)) 



There is also a truth function which gives the current truth value 
associated with a TMS node. Thus (truth (mammal 'Joe)) would evaluate to true 
after the above assertion had been made. In other circumstances it might 
evaluate to false, or unknown. There is also a why primitive which takes a TMS 
node and gives an explanation for its truth value in terms of other truth values of 
TMS nodes which imply that value. 

The use of TMS predicates as a knowledge representation mechanism 
results in comprehensible modular constructs which interface cleanly to a truth 
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maintenance system. 



Canonical Naming 

The central theme of this document is equality. A TMS predicate has 
been defined, == which takes two designators and returns a TMS node 
representing the fact that the two designators have the same referent. As was 
mentioned earlier my use of the word "designator" is equivalent to the use of the 
word "term" in systems of formal logic. The word "designator" is used here for 
no better reason than that it focuses on the distinction between a term and its 
referent. The == predicate interprets lisp s-expressions as designators and 
designators will be uniformly represented by lisp s-expressions throughout. 

One basic problem with the existence of more than one designator for an 
individual object is that knowledge about that individual is often only given in 
terms of a single designator. For example consider the situation defined by the 
following assertions: 

(assert (== ' (residence the-US-president) 
' the-whi te-house) ) 

(assert (== 'Jimmy-Carter 'the-US-president)) 

Now suppose we ask for the truth of an equality between (residence 
Jimmy-Carter) and the-white-house. We would like the system to see that since 
Jimmy Carter is the U.S. president, his residence is the white house. There must 
be some way of merging the knowledge about a single object or individual which 
is given in terms of the various designators for that object or individual. This 
merging of knowledge can be accomplished via a canonical name for each 
individual. Any knowledge about an object can be stated in terms of its 
canonical name by substituting canonical names into expressions. This allows 
knowledge to interact with other knowledge which is similarly stated in terms of 
the canonical name. 

The canonical name strategy for handling equality has been used in 
many previous problem solvers. The primary innovation here is the use of this 
mechanism as a primary reasoning strategy as will be discussed in the next few 
sections. In the present system only knowledge in the form of equalities fully 
benefits from canonical naming since substitutions are only performed on 
designators and not on TMS assertions. The details of the substitution process 
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will be given in a later section; the maintenance of canonical names being 
concentrated on here. 

The == predicate creates auxiliary data structures which are used in 
maintaining canonical names for objects. A set of equalities defines an 
equivalence relation on designators. Since all of the designators in an equivalence 
class have the same referent, one of those designators can be chosen as the 
canonical name for that referent. Only equalities whose TMS nodes are true are 
used in this process and the canonical names are incrementally maintained when 
the truth values associated with equalities change. There is a "what-is" function 
which takes a designator and returns the canonical name for that designator to 
the user. This function is surprisingly useful, as will be seen in the next few 
sections. 

Algebraic Unknowns 

Designators, and the concepts involved in equality are related to the 
concepts involved in the use of algebraic unknowns. Notice that algebraic 
unknowns can be separated from bound variables and arbitrary individuals which 
can be involved in instantiation processes. The ways in which algebraic 
unknowns take on values is entirely independent of such instantiation processes. 
Consider a system of n independent linear equations in m unknowns. If there are 
more unknowns than there are equations, then there are several degrees of 
freedom left in the system. However, it is possible to "eliminate" n of the 
unknowns by expressing each of them as a linear combination of the others. The 
"value" of each of these eliminated unknowns could then be set to the 
corresponding linear combination of the remaining unknowns. 

It is possible to separate the issue of unknown elimination from algebraic 
manipulation. Algebraic manipulation is the deduction of new equalities from 
given equalities. For example, from the equality y = (+ (* a x) b) it is possible 
to deduce (* a x) = (- y b). Neither of these equations need have anything to 
do with the elimination of unknowns; they are simply assertions of equality. I 
will attempt here to relate the notion of an algebraic unknown, the act of giving 
unknowns values, and the use of equality, via a general notion of canonical 
naming. 

Under the view of unknowns presented here, they are a type of 
designator. I am not immediately concerned with which designators qualify as 
unknowns but will deal only with the general notion of a designator. In light of 
the above discussion of canonical names the act of giving an unknown a value 
can be viewed as the act of choosing a designator, other than that unknown, as 
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the canonical name for the unknown's class. For example if y and (+ (* a x) b) 
are in the same class, i.e. y = (+ (* a x) b), and for some reason (+ (* a x) b) 
is chosen as a better designator than y to represent the referent of the class, then 
y can be said to have the value (+ (* ax) b). 

The next few sections will develop propagation of constraints as an 
algorithm for solving systems of equations. It is shown that this algorithm can be 
viewed as a process of shifting canonical names. 

Constraint Propagation as Computation via Naming. 

Lately there has been much interest in, and development of, constraint 
propagation as an efficient algorithm for reasoning about sets of mutually 
constrained quantities [Steele & Sussman 78]. This algorithm is also interesting in 
that it seems to simulate one of the ways human engineers reason. In general a 
constraint propagation system has a set of "cells" which can take on values, and a 
set of "constraints" which constrain those values. In constraint propagation 
whenever a deduction can be made from the previously determined values and a 
single constraint this deduction is made. Each such deduction assigns a new 
value to a cell. In what will be termed "simple" constraint propagation these are 
the only deductions which are made. Constraint propagation terminates when 
there are no deductions which can be made from the cell values and a single 
constraint. It is easy to see that such a deduction process can take no longer 
than linear time in the number of constraints. 

A method of implementing constraint propagation is presented here 
which is based entirely on the algorithms for handling equality discussed above. 
The "cells" of traditional constraint propagation are implemented as equivalence 
classes of designators which take on "values" in the form of canonical names for 
those classes. Each constraint is implemented as a set of equalities, and the 
primary deductive mechanism employed is substitution controlled by a method of 
choosing the canonical names of equivalence classes. The constraint propagation 
techniques developed here are relied upon heavily in an electrical circuit analysis 
system described in a later chapter. 

Constraints are created via constraint predicates. The basic constraint 
predicates are +constrained and *constrained, which are defined as follows: 
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(defpred +constrained (sum al a2) 
(== sum * (+ ,al ,a2)) 
(== al * (- .sum ,a2)) 
(== a2 '(- .sum .al))) 

(defpred *constrained (prod ml m2) 
(== prod ' (* .ml ,m2)) 
(-> (prednot (== m2 0)) 

(== ml '(// , prod , m 2))) 
(-> (prednot (== ml 0)) 

(== m2 '(// ,prod ,ml))) 
(presumably (prednot (=» ml 0))) 
(presumably (prednot (=- m2 0)))) 

The backquote macro is used to simplify the creation of designator s- 
expressions. The backquote macro is a form of quote in which items in the 
interior of the backquoted structure preceded by a comma are replaced by their 
value. Thus '(+ ,x ,y) is equivalent to (list '+ x y). 

Note the use of the "presumably" predicate. If a *constrained is in 
force (its TMS node is true) then the node representing (== ml 0) will default to 
false. However if this leads to a contradiction (ml is discovered to be 0), then 
the support for ml not being given via the default construct can be 
automatically retracted by the TMS. 

Using the constraint predicates defined above, it is possible to construct 
constraint networks within the equality system. Given an appropriate method of 
deciding which designator of two given designators is "better" for use as a 
canonical name, the constraint nets formed in this way can interact with the 
canonical equality system to actually propagate values. To decide when one 
designator is a better designator than another it is necessary to define a partial 
order on designators which will be referred to here as the better-name relation. 
Two designators can be unordered by the better-name relation and in that case 
whichever becomes the canonical name first will remain so, thus avoiding undue 
substitution. 

In this discussion of constraint propagation a better-name relation is 
assumed in which designators are divided into two classes, "known" and 
"unknown". For atomic designators a predicate is provided to determine whether 
the designator is known or unknown, and a functional expression is known iff its 
operator and all of its arguments are known. Known designators are always 
better than unknown designators. This convention along with substitution is 
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enough to give propagation of constraints. 

The two methods by which designators are generated automatically in 
the equality system are substitution and evaluation. Only substitution will be 
considered here; evaluation (which is the way algebraic simplification gets done) 
will be discussed later. Whenever a new canonical name is assigned to a class 
that canonical name is substituted into designators which contain other members 
of that class. Such substitutions generate new designators from a pre-existing 
ones. Whenever this is done an equality is created between the original 
designator and the one resulting from the substitution. The truth of this equality 
is implied by the truth of the equalities used in the substitution. A more detailed 
description of the implementation of this process will be given in chapter III. 
Now consider the following constraint net taken from [Steele & Sussman 78]: 



skfc> 



pi 



Figure 2. A Simple Constraint Network. 

This constraint net can be created in the equality system by the following 
assertions: 

(assert (*constrained 'pi 'ml 'm2)) 

(assert (+constrained 'si 'al *a2)) 

(assert < — 'pi 'al)) 

(assert (== *a2 'ml)) 

The initial equivalence classes are derived from the equalities and a 
canonical name is chosen for each class. I will assume that all designators used 
so far are classified as unknown and that, all other things being equal, smaller 
designators are a better choice for a canonical name than larger, more complex 
ones. The equivalence classes defined so far and their (somewhat arbitrary) 
canonical names are: 



Chapter I 13 Canonical Naming 



m2, (// pi ml) m2 

ml, (// pi m2), a2, (- si al) ml 

pi, (* ml m2), al, (- si a2) pi 

si, (+ al a2) si 

The substitution process would add several equalities which are shown 
below (each equality has been associated with the canonical name for the class 
which contains its arguments). 

(== ' (+ al a2) *(+ pi ml)) si 

(== ' (- si al) '(- si pi)) ml 

(== ' (- si a2) M- si ml)) pi 

Since none of the generated designators are known, the original 
canonical naming is stable and no propagation occurs. Now suppose however 
that si is made equal to 4 and al is made equal to 2. 4 and 2 are assumed to be 
known and therefore the canonical designators of the classes corresponding to si 
and al are changed to 4 and 2 respectively. Substitution then generates the 
following equalities: 

= ' (+ pi ml) '(+ 2 ml)) 4 

= M- si ml) '(- 4 ml)) 2 

= '(- si pi) '(- 4 2)) ml 

= '(// pi ml) '(// 2 ml)) m2 

= ' (// pi m2) '(// 2 m2)) ml 

As the substitution process was described above, other substitutions 
would also be made (into (+ al a2) for example). However these designators 
have already been substituted into. The designator generated by the earlier 
substitution is equal in all its parts to the original designator, and can be further 
substituted into, making substitution into the original designator redundant. The 
details of the way in which the system avoids redundant substitutions will be 
given in chapter IV. 

Given the above equalities a2 has been made equal to (- 4 2) since it 
was originally equal to (- si al). Since (- 4 2) is a known designator it must 
become the canonical name of a2's class. This change in canonical name then 
leads to the following further equalities via substitution: 
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(== ' (+ 2 ml) M+ 2 (-4 2))) 4 

(== • (_ 4 ml) ' (- 4 (-4 2))) 2 

(== ' (* ml m2) ' (* (-4 2) n>2)) 2 

(== ' (// pi ml) M// 2 (- 4 2))) m2 

Now notice that m2 has been made equal to (// 2 (- 4 2)). Since the 
new value is known it becomes the canonical name and further substitution leads 
to no changes in canonical names. Let the "value" relation x=>y be defined as 
meaning y is the canonical name of the equivalence class of x. The results of the 
above propagation can be summarized by the following relations: 



al 


=> 


2 


Pi 


=> 


2 


si 


=> 


4 


a2 


=> 


(- 4 2) 


ml 


=> 


(- 4 2) 


m2 


=> 


(// 2 (- 4 2)) 



If m2 and si had been specified instead of al and si, no propagation 
would have occurred even though all values are completely determined. Such 
situations are discussed by Steele and Sussman and are due to loops in the 
constraint net. One technique which has been successfully used to handle 
involves propagating symbolic values which are classified as known, just as 
numeric values are propagated. This technique, called plunking, is discussed later 
in its own subsection. 

A general feeling for the efficiency of this algorithm can be gotten by 
considering the number of substitutions performed. As mentioned above once a 
new designator has been derived via substitution the designator which was 
substituted into will not be substituted into again (assuming no equalities are 
retracted which invalidate the original substitution). Therefore each designator 
gives rise to a string of descendants generated by substitution, each descendant 
generating the next. The number of such descendants for a given designator can 
be no larger than the sum of the number of times the top level subdesignators 
change canonical names. The number of top level subdesignators is one plus the 
arity of the function involved, and is assumed quite small. The number of times 
the canonical name of a class changes can also be assumed to be quite small in 
practice, since it is usually due to a shift from an unknown to a known 
designator. Thus the total number of substitutions is bounded by a small 
constant times the number of original designators in the constraint net. The 
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evaluation of designators, which has not been discussed yet, can generate 
designators in ways not involving substitution. These other designators expand 
the set of names which are substituted into and can make constraint propagation 
somewhat more expensive. 

It should be emphasized that the propagation process is really driven by 
the decisions about which designators are better canonical names. In chapter III 
the discussion of electronic circuit analysis will require a more elaborate better 
name relation which relies on a destinction between referentially transparent and 
referentially opaque designators. 
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Chapter II 

Procedural Embedding 

The three sections of this chapter deal with various ways in which 
procedures have been embedded in the equality system and the ways in which 
they are used. The first section introduces the major mechanism used in 
procedural embedding and shows how simple integer arithmetic can be done using 
this mechanism. The next section discusses some problems surrounding equalities 
which should have false truth values, such as equalities between distinct integers. 
The final section deals with the problem of solving for values which are 
determined by a set of constraints but can not be solved for directly via 
propagation of constraints. 

Evaluation Functions 

Often quantified knowledge can be associated with a specific function or 
operator. Knowledge about the properties of algebraic operators is an example of 
this type of knowledge. To facilitate the control of the instantiation of such 
knowledge a special mechanism involving "evaluation functions" has been 
incorporated into the equality system. This section will first develop the general 
mechanisms involved in evaluation functions and then turn to applications 
involving algebraic operators. 

Consider the knowledge that the mother of every animal is female. This 
information is associated with the mother function since it only applies to 
applications of that function. Formally such knowledge would be represented as 
Vx(-> (animal x) (female (mother x))) which could then be instantiated with any 
term. Knowledge concerning a certain operator can be attached to that operator 
in the equality system in the form of evaluation functions. Each evaluation 
function for a given operator can be applied to any designator representing an 
application of that operator and performs an instantiation of the quantified 
knowledge it contains. The above knowledge about the mother function could be 
embedded in the following evaluation function: 



Chapter II 17 Procedural Embedding 



(def pred animal (x) ) 
(defpred female (x) ) 

(defun mother-eval-1 (designator-s-expression) 
(if (eq (car designator-s-expression) 'mother) 
(let ((x (cadr designator-s-expression))) 

(assert (-> (animal x) (female designator-s-expression) )))) ) 

This function must be appropriately associated with the mother operator 
in the equality system for it to be useful. To simplify things a special mechanism 
has been created to create such functions and associate them with operators. 
Using this mechanism the above function could have been both defined and 
associated with the mother operator as follows: 

(eval fun mother (x) 

(-> (animal x) (female '(mother ,x)))) 

As in the case of a defpred, the body of a evalfun is a list of forms 
which return TMS nodes when evaluated. The function created by a evalfun 
evaluates the forms of the body in an environment in which the argument list of 
the evalfun has been bound to the appropriate parts of a designator representing 
an application of the involved operator. The nodes created by the forms of the 
body are asserted as premises. 

There can be more than one evaluation function for a given operator. 
For example the knowledge that the mother of an animal is also a parent of that 
animal could be represented in an evaluation function as: 

(defpred parent (x y)) 

(eval fun mother (x) 

(parent '(mother ,x) x)) 

The control of the application of evaluation functions (and therefore the 
control of the instantiation of the quantified knowledge they contain) will be 
discussed in detail in chapter four. For now however there is no problem in 
assuming that it is simple reference instantiation, that is to say that the 
quantified knowledge is instantiated (i.e. the evaluation functions are applied) 
whenever reference is made to some application of the involved operator. 

This mechanism can be used to achieve simplification of numerical 
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expressions in the equality system. Consider the following: 

(defun all-numbers (things) 
(or (nul I things) 

(and (numberp (car things)) 

(all-numbers (cdr things))))) 

(evalfun + addends 

(if (a I I -numbers addends) 

(== ' (+ , ©addends) (apply 'plus addends)))) 

(eval fun - args 

(if (a II -numbers addends) 

(-- '(- .eargs) '(apply 'minus args)))) 

etc. 

Here a slightly different syntax has been used to handle a variable 
number of arguments in which the argument list has been replaced by a single 
atom which is bound to a list of the designators representing the arguments of 
the operator. Also the form in the body of the evalfun can return nil, in which 
case nothing is done (this is also true in the body of a defpred). The use of ",@" 
m the interior of a backquote explodes the value of the form following (which 
must be a list) up into the list in which the ",@" appears. For example '(a ,(list 
'a »b) ,@(list 'a 'b)) evaluates to (a (a b) a b). 

Using these definitions, evaluation of the designator (- 4 2) would result 
in an equality being added between it and 2. It is up to the better-name relation 
to determine that 2 is a better designator for the canonical name of the 
equivalence class than (- 4 2). If 2 becomes the canonical name of its class, then 
it will be substituted into other designators such as (// 2 (- 4 2)). The results of 
such substitutions can then be evaluated to yield further simplifications. The 
definition of an evaluation function for // must be careful not to state false 
equalities when two integers do not divide evenly. 

It should be clear from the way the LISP plus and minus functions were 
used above that arbitrary procedures can readily be embedded into the evaluation 
functions. Special purpose algebraic simplification routines are embedded in the 
equality system in exactly this way. Thus an evaluation functions for +, -, etc. 
have been created which are useful when the arguments to the operator are 
symbolic. These evaluation functions create an equality between the expression 
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being evaluated and a symbolically simplified (or standardized) form of it. The 
specific nature of the algebraic simplification algorithms used is discussed in 
appendix two. 

Disequalities 

As was mentioned in the section on constraint propagation, it is possible 
that a set of constraints determine values for designators, but that constraint 
propagation does not result in those values being found. This situation results 
when loops are present in the constraint net and is illustrated by the following 
example: 

(♦constrained 'y 2 ' x) fy = 2x1 
(+constrained 3 'x *y) ix+y = 3} 

Notice that even though the value of x and y are determined by the 
constraints, no constraint propagation occurs since all of the designators for x and 
y contain either x or y and are therefore classified as unknown (assuming x and y 
are classified as unknown). How the equality system actually finds values which 
are determined but not deduced via constraint propagation is the subject of the 
next section. This section is only concerned with the possible contradictions that 
can arise in these situations. 

Suppose that an equality between x and 3 is added to the above 
constraints. Propagation would occur and the system would deduce that y equals 
(* 2 3) by substituting into (* 2 x). It would also deduce that y equals (- 3 3) 
by substituting into (- 3 x). The evaluation functions attached to the operators 
would then put both and 6 into the equivalence class for y. Nothing in the 
system so far says that this is impossible (the names being used for numbers 
might be representing numbers in modular arithmetic, and indeed does equal 6 
if they refer to equivalence classes mod 3). 

To rule out such interpretations of numeric designators in the equality 
system a special check for numeric designators has been placed in the == 
predicate. Where the designators passed to == are distinct numbers == insures 
that the node it returns to represent the equality is false. The equality system 
internally creates equality nodes between every designator in an equivalence class 
and the canonical name of that class. Thus if a numerical designator is the 
canonical name for some class the entrance of any other numerical designator 
into that class will result in a contradiction via the equality node which would be 
created between the two numerical designators. 
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This technique could easily be generalized to a method for ensuring the 
integrity of other disequalities. Certain designators could be called "truely 
canonical" designators in that a "truely canonical" designator must always be the 
canonical name for its equivalence class. Thus no two "truely canonical" 
designators could ever be equal since then both of them would have to be the 
canonical name of the same class. Therefore whenever an equality assertion is 
created between two "truely canonical" designators the system could force that 
equality to be false. The above method for handling numbers would then be 
implemented by simply making numeric designators "truely canonical", but the 
method is certainly not restricted to numbers and could be used to ensure that, 
for example, Earth is not the same planet as Mars. 

Plunks 

One method for finding values which are determined but not found via 
constraint propogation is referred to as "plunking" and is discussed by Stallman 
and Sussman [Stallman & Sussman 77]. In their system plunking involves giving a 
variable a symbolic value which propagates through the network just as any 
numeric value would. Certain equalities generated during this propagation can 
then be passed to an equation solver which solves for the plunks in terms of 
other quantities. 

The algorithm discussed here is slightly different from the one used in 
their constraint language due to the difference in the representation of the 
constraints. In the equality system a class which contains only unknown 
designators is chosen to be plunked. A new unique designator is generated to 
represent a symbolic value, which will be called a plunk. This designator is then 
made equal to some designator in the plunked class. For the plunking to work 
properly the better name relation must treat the plunks specially. A plunk is 
treated as a known designator, and it induces propagation. However designators 
containing plunks are considered worse names than other known designators so 
that designators with plunks are replaced by designators without them if the 
opportunity arises. The better name relation as it has been described so far 
divides designators into three classes and orders those classes as follows: 

known designator without plunks 
are better than: known designators with plunks 

are better than: unknown designators 
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Now consider the constraint net described in the previous section: 

y = 2x 
x+y = 3 

In this constraint net it is possible to plunk x by adding an equality 
between x and a designator generated to represent the plunk, call it plunk- 1. The 
canonical name for x's equivalence class would shift to plunk-1 which would in 
turn induce substitutions, y is equal to (* 2 x) which generates (* 2 plunk-1) 
under the substitution. Since this later designator is treated as known the 
canonical name for y's class becomes (* 2 plunk-1). But y is also equal to (+ x 
1) which generates (+ plunk-1 1) when substituted into. 

To solve for the plunks the notion of a coincidence has been adapted 
from previous constraint systems to use in the equality system. A coincidence 
occurs when a "known" designator enters a class whose canonical name was 
already a "known" designator. In such cases the pair of known designators is 
passed to a coincidence handler. The coincidence handler checks to see if they 
are simply substitution variants of each other and if not it then checks to see if 
either contain a plunk. If a plunk is present then the coincidence handler 
attempts to use the coincidence to solve for a plunk. In the above example the 
coincidence handler would be passed (* 2 plunk-1) and (+ plunk-1 1). Since the 
pair of designators which are passed to the coincidence handler are equal it can 
solve this equality for plunk-1 using a solver from the algebraic simplification 
system. The coincidence would now add an equality between plunk-1 and 1. 
Since known designators which do not contain plunks are better than known 
designators containing plunks, 1 would become the canonical name for plunk- l's 
equivalence class, and 1 would be substituted for plunk-1 wherever it occurred. 

In general more than one plunk may be necessary to solve a set of 
constraints. In that case the designators which are passed to the coincidence 
handler may contain more than one plunk. To handle such situations each plunk 
is given a plunk-weight which is used by the better-name relation. A lower 
weight plunk is always considered a better name than a higher weight plunk, and 
the plunk-weight of a designator containing no plunks is considered to be zero. 
The plunk-weight of designators containing plunks is the maximum weight of the 
plunks in the designator. Now when the coincidence handler gets expressions 
which contain more than one plunk it solves for the highest weight plunk in 
terms of the others. Since any designator containing only the other plunks will 
have a lower plunk weight, the canonical name of the solved for plunk must shift, 
and therefore that plunk will be replaced wherever it occurs. 
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Plunking is not performed automatically in the present equality system. 
To invoke plunking a soive-for function has been defined which takes a single 
designator and plunks quantities which that designator can be expressed in terms 
of. The net effect is that the argument passed to solve-for is made equal to a 
"known" expression which may contain plunks. If plunks are present in the 
resulting expression the solve-for function can be used again on the plunks 
themselves. This can be done iteratively to solve for in expression which is 
indeed determined by a set of constraints. The following dialogue shows how 
this procedure can be applied to the above constraint net. 

(assert (*constrained 'y 2 'x)) 
(♦CONSTRAINED 'Y 2. *X) 

(Assert (+constrained 3 'x 'y)) 
(+CONSTRAINED 3. 'X *Y) 

(uhat-is *x) 
X 

(uhat-is 'y) 
Y 

(solve-for 'x) 
1 

(uhat-is 'y) 

2 ' 

The next chapter applies the mechanisms which have been developed 
so far to the domain of electrical circuit analysis. 
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Chapter III 

Application to Electronics 

In this section a system for analyzing electronic circuits is presented 
which has been implemented using the equality system. Only a minimum of 
knowledge about electronics is necessary for the comprehension of this 
section and the predicates defined here contain all the electronics knowledge 
used. Since many general properties of the equality system are exhibited, it 
is hoped that even readers with no background in electronics will take the 
time to read this section. 

A set of basic electronic predicates are defined, which can be used 
to build arbitrary circuits. However, for the constraint propagation to work 
properly in reasoning about circuits defined with these predicates, a better- 
name relation must be used which is slightly more complex than those 
discussed so far. A distinction is made between referentially transparent and 
referentially opaque designators, the former being a better name than the 
latter. The need for this distinction, and its use, is discussed in the context 
of circuit analysis. 

Basic Predicates 

The predicates which are defined below are used to create some of 
the constraints associated with various device types which the electrical 
analysis system deals with. Predicates are also given which are used to wire 
the device components together. The predicate c== has been written in lisp 
as a special constraint predicate. c== takes two designators and solves the 
implied equation for each internal designator in the two expressions. This 
results in a set of equalities which are all implied by a node returned by c==. 
Appropriate quantities are assumed not to be equal to zero as was done in 
♦constrained. (c== 'a '(* b c)) is equivalent to (*constrained 'a 'b 'c) but 
less efficient. 

The use of » in designators gives a shorthand for the repeated 
functional composition of monadic functions. Thus (» voltage rl circuit3) 
is treated identically to (voltage (rl circuit3)). This is often very convenient 
when dealing with complex structured objects. 

(defpred inv-constrained (a b) 
(== a '(- ,b)) 
<== b '(- ,a))) 
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(defpred one-port (element) 

(+constrained ' (» potential tl .element) 
' (>> potential t2 .element) 
' (voltage .element)) 
(== '(current .element) 

'(>> current tl .element)) 
(inv-constrained '(» current tl .element) '(current t2 .element))) 

(defpred resistor (r) 
(one-port r) 
(♦constrained '(voltage ,r) '(current ,r) '(resistance ,r))) 

(defpred voltage-source (vs) 
(one-port vs) 
(== '(voltage ,vs) '(strength ,vs))) 

(defpred current-source (cs) 
(one-port cs) 
(== '(current ,cs) '(strength ,cs))) 

There are no definitions for controlled sources given here. This is 
because a controlled source can be simply represented as a one-port with an 
appropriate current or voltage constraint. 

The most general way to handle wiring is to write a "node" predicate 
which takes any number of terminals and creates equalities stating that that the 
potentials of all the terminals are the same and that the sum of the currents is 0. 
This predicate definition uses a slightly different syntax since it is a predicate on 
an arbitrary number of arguments. The atom representing its argument list is 
bound to a list of the values of the arguments in a call to the predicate. The 
two other wiring predicates given are a little more natural to use than the node 
predicate. 

(defpred node terminals 

(c== ' (+ ,@(mapcar '(lambda (term) '(current .term)) terms))) 
(let ((node-potential '(potential .(car terms)))) 
(predand (mapcar '(lambda (term2) 

(=« node-potential '(potential .terml))) 
(cdr terms))))) 
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(defpred connected (tl t2) 

(== '(potential ,tl) '(potential ,t2)) 
(== '(current ,tl) M- (current ,t2)))) 

(defpred exists (x) ) 

(evalfun composite (tl t2) 

(-> (exists '(composite , tl ,t2)) 

(predand (== '(potential ,tl) '(potential ,t2)) 

(== '(potential (composite , tl , t2) ) '(potential ,tl)) 
(c== '(current (composite ,tl ,t2)) 

'(+ (current ,tl) (current , t2))))) 
(presumably (exists '(composite , tl ,t2)))) 

These predicates can be used to create circuits. As a trivial first 
example a circuit which deals only with a single voltage source and resistor is 
defined below and shown in figure 3. 

(defpred ohm-test (circuit) 

(voltage-source *(vs .circuit)) 

(resistor ' (r .circuit)) 

(node ' (» tl vs .circuit) ' (» tl r .circuit)) 

(node '(» t2 vs .circuit) '(» t2 r .circuit))) 




Figure 3. The "Ohm-test" Circuit. 
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A particular instance of this circuit is created by asserting an application 
of the above predicate to a specific circuit name. The following scenario shows a 
simple use of this predicate. 



(assert (ohm-test 'cl)) 
(OHM-TEST Cl) 

(why (one-port ' (r CD)) 
((ONE-PORT (R CD) IS TRUE FROM 
(1 (RESISTOR (R CD) IS. TRUE)) 

(why (resistor ' (r cD ) ) 
((ONE-RESISTOR (R CD) IS TRUE FROM 
(1 (OHM-TEST CD IS TRUE)) 



This method of circuit definition, while quite clean and semantically 
satisfying has the disadvantage that the circuit topology of an instance of the 
circuit definition can not be incrementally altered because the assertion that the 
circuit is an ohm-test completely determines the structure of the circuit. Such 
incremental alterations would be possible if the various assertions about the 
circuit were taken as independent premises as was done by Stallman and Sussman 
in their electronic analysis system ARS [Stallman & Sussman 77]. This can be 
done by using a simple lisp function to construct the circuit which might be 
defined as follows: 

(defun make-ohm-test (circuit) 

(assert (voltage-source ' (vs .circuit))) 

(assert (resistor Mr .circuit))) 

(assert (node '(» tl vs .circuit) * (» tl r .circuit))) 

(assert (node '(» t2 vs .circuit) '(» t2 r .circuit))) 

t) 
MAKE-OHM-TEST 

(make-ohm-test 'cD 
T 

(why (restistor Mr cl))) 

((RESISTOR (R CD) IS TRUE AS A PREMISE) 
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Now it is possible to change the topology of this circuit by incrementally 
altering the basic premises which determining that topology. For example the 
following scenario shows how another resistor could be spliced into the circuit: 

(assert (resistor '(r2 cl))) 
(RESISTOR R2) 

(retract (node '(» tl vs cl) '(tl r cl))) 
T 

(assert (node '(» tl vs cl) '(» tl r2 cl))) 
(NODE (» Tl VS Cl) (» Tl R2 Cl)) 

(assert (node '(» t2 r2 cl) '(» tl r cl))) 
(NODE (» T2 R2 Cl) (» Tl R Cl)) 

In the discussions of electronics that follows, circuits will be defined via 
the defpred mechanism but the reader should be aware that circuits which can be 
incrementally modified can be constructed and that derived knowledge about 
such circuits is also incrementally modified when changes are made. 

The Better-Name Relation in Electronics 

The primary reasoning strategy used in the electronic analysis system 
described here is propagation of constraints. As propagation of constraints was 
defined above it is driven by the better-name relation, that is the ordering on 
designators which determines when one designator is preferred over another as 
the canonical representative for an equivalence class of designators. Here an 
investigation is made into the properties a better-name relation should have such 
that constraint propagation is done in a useful manner. 

Consider the simple circuit defined in the previous section and shown in 
figure 3. A trivial problem in circuit analysis would be to make an ohm-test 
circuit, call it cl, and ask for the current in the resistor. The current in the 
resistor has many designators, the simplest of which is (» current r cl). 
However the designator which we would like to get as a value for the current is 
the ratio of the strength of the voltage source to the resistance of the resistor. 
The better name relation should be such that designators in terms of source 
strengths and component parameters are better names that those in terms of 
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currents and voltages. 

An interesting observation is that designators involving voltages and 
currents are sensitive to the context in which they occur. For example the 
voltage across a certain resistor is not purely a function of that resistor but is 
sensitive to the environment of the resistor. A designator whose referent is 
context sensitive is termed referentially opaque, otherwise it is termed 
referentially transparent. It seems that referentially transparent designators are 
better names, at least in the context of electronic analysis. The equality system 
presently perfoms a simple syntactic search for potential, current, or voltage 
to determine if a designator is opaque or transparent. 

Numeric designators should be better names than transparent but non- 
numeric ones. Thus a numeric value for a component parameter is a better-name 
than a designator in terms of the parameter function applied to that component. 
If plunks are present then they are considered to be transparent so that 
propagation treats them properly. The classification of designators by the better 
name relation is now as follows: 

numerical designators 
are better than: transparent designators uithout plunks 
are better than: transparent designators nith plunks 
are better than: opaque designators 

The better-name relation for the electronics system also takes the 
designator size into account when the designators are of the same type, smaller 
designators being better. Given this framework it is now possible to see how the 
equality system analyzes circuits. 

A Simple Analysis 

The following simple dialogue represents an analysis of an instance of 
the simple ohm-test circuit defined above. In all the sample dialogs presented 
here things typed by the user appear in lower case while responses typed by the 
system appear in upper case. 

(assert (ohm-test 'cl)) 
T 

(what-is ' (>> current r cl)) 

(// (» STRENGTH VS Cl) (» RESISTANCE R CD) 
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To see how this analysis took place it is necessary to examine the 
original equalities which were provided to the system. Among these equalities are 
the following four: 

(== ' (>> potential tl r cl) 

' (>> potential tl vs cl)) 

(== ' (>> potential t2 r cl) 

' (>> potential t2 vs cl)) 
(== ' (>> voltage vs cl) 

' (- (potential tl vs cl) 

(potential t2 vs cl))) 

(== * (>> vol tage r cl) 

* (- (potential tl r cl) 

(potential t2 r cl))) 

The first two equalities are created by the node predicate. A canonical 
name will be chosen for the potential of each node. Since the alternatives for 
these canonical names are roughly equivalent in their type and complexity, the 
choice is arbitrary, but let us assume that the canonical names are given in terms 
of the voltage source. So (» potential tl vs cl) is the canonical name for (» 
potential tl r cl) and similarly for the other potential. This choice of canonical 
names results in the following equality being generated via substitution: 

(== ' (- (» potential tl r cl) (» potential t2 r cl)) 
' (- (» potential tl vs cl) (» potential t2 vs cl))) 

This equality interacts with the last two of the first four equalities above 
making the voltage of the voltage source equal to the voltage of the resistor. By 
the definition of a voltage source, the voltage of the source is equal to the 
strength of the source which has a referentially transparent designator. The 
designator (» strength vs cl) therefore becomes the canonical name for the 
voltage of both the voltage source and the resistor. Now consider a result of 
Ohm's law which was generated as a part of the constraints in the definition of a 
resistor: 

(== ' (>> current r cl) 

' (// (voltage r cl) (resistance r cl))) 
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Substitution will yield the equality: 

(== ' (// (voltage r cl) (resistance r cl)) 
' (// (strength vs cl) (resistance r cl))) 

The ratio of the strength of vs to the resistance of r is a transparent 
designator, and therefore becomes the canonical name for the current in the 
resistor and the analysis is complete. Further discussions about circuit analysis 
will be a little less formal about the actions of the equality system. Instead of 
referring to the addition of a referentially transparent designator to an 
equivalence class, it will simply be said that designators in that class have been 
determined. Thus, instead of saying that a referentially transparent canonical 
name for the current in r has been found, I will simply say that the current in r 
has been determined. 

Unconstrained Plunks 

There are many circuits which cannot be analyzed using the definitions 
given so far without resorting to plunks. In many circuits the need for plunks is 
quite surprising since human engineers can easily analyze the circuit by 
inspection. The last three sections in this chapter discuss ways in which the 
electrical analysis system can be made more powerful. In this section the notion 
of an "unconstrained plunk" is developed for use in situations in which there are 
unconstrained degrees of freedom left in the quantities being reasoned about. 
The next section discusses some of the ways in which multiple descriptions can be 
used to aid reasoning in constraint propagation. The final section discusses the 
effect of restructuring constraints in ways which correspond to "redrawing" 
circuits. 

As a first example of a circuit in which constraint propagation alone 
fails to perform the desired analysis, consider the circuit shown in figure 4. 
When an instance of this circuit is analyzed the current in each resistor is 
determined by the current through the current source. The current in each 
resistor determines the voltage across that resistor so it would seem that the 
voltage across the current source should also be determined by the system. This 
however does not happen since the voltage across the source is equal to the 
difference between the potentials of its two terminals, and neither of these 
potentials can be determined. Nor is the difference determined directly since 
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Figure 4. The "Series-1" Circuit. 

(defpred series-1 (circuit) 

(current-source * (cs .circuit)) 

(resistor ' (rl .circuit)) 

(resistor ' (r2 .circuit)) 

(series-1-topology ' (cs .circuit) 4 (rl .circuit) * (r2 .circuit))) 

(defpred ser ies-1-topology (cs rl r2) 
(node *(tl ,cs) 4 (tl ,rl)) 
(node Mt2 ,rl) Mtl ,r2)) 
(node Mt2 ,r2) Mt2 .cs))) 

there is no statement that the voltage across the source is the sum of the voltages 
across the resistors. If there were some way in which potentials could be 
determined, instead of potentials differences, this problem could be solved. 
However there is a degree of freedom left undetermined in the node potentials 
and none of these potentials can be solved for. 

In the framework of circuit analysis presented so far the above circuit 
could be ^analyzed via plunking. In such a scheme some node potential could be 
"plunked" by setting it equal to a plunk designator. Propagation would then take 
place and each node potential would take on a canonical name in terms of the 
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plunk. The voltage across the source would then be determined in terms of the 
plunk, which would cancel out of the resulting expression. Previous electrical 
analysis programs which relied on propagation of constraints have used plunking 
in this way [Stallman & Sussman 77] [de Kleer & Sussman 78]. 

In the use of a plunk to analyze the above circuit the plunk cannot be 
eliminated since the absolute values of the node potentials are not determined. 
Plunks which can never be eliminated will be called "unconstrained plunks". The 
special treatment the system gives to ordinary plunks by the better-name relation 
and the coincidence handler is designed to ensure that plunks can be eliminated. 
Since an unconstrained plunk cannot be eliminated, it need not be given this 
special treatment. In fact no distinction need be made whatsoever between this 
type of plunk and ordinary referentially transparent designators. However, since 
no attempt is made by the system to solve for unconstrained plunks, it is 
important that such unconstrained plunks be mutually independent, that is to say 
that no one of them is expressible in terms of the others. 

Unfortunately the fact that a plunk can never be eliminated is an 
observation made outside the system (the equality system is not capable of 
realizing this simply from a given set of equalities). However since an 
unconstrained plunk can be treated exactly as a simple transparent designator the 
"plunk" can be made explicitly by the user of the equality system simply by 
stating an equality between the unconstrained quantity to be plunked and a 
referentially transparent designator which will act as the plunk. This gives the 
user more control over the reasoning process by explicitly stating a quantity to be 
plunked. The choice of which quantity is given an unconstrained plunk in a 
system of mutually constrained quantities can also be important to the human 
user who wishes to see results in terms of certain quantities. 

In all electronic circuits there is at least one degree of freedom in the 
node potentials which is not determined by the circuit. This degree of freedom 
can be dealt with by defining a reference node to be given an unconstrained 
plunk. For convenience the node given the unconstrained plunk will be called 
the ground node. The ground is specified simply by setting its potential equal to 
the designator "ground-potential" (the unconstrained plunk) which is treated as a 
referentially transparent designator. To see how the specification of a ground 
node interacts with the analysis of the above circuit consider the following 
alternate definition of the circuit's topology: 
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(defpred ser ies-1-topology (cs rl r2) 
(node '(tl ,cs) Mtl ,rl)) 
(node Mt2 ,rl) Mtl ,r2)) 

(== M>> potential t2 ,r2) 'ground-potential) 
(== M» potential t2 ,cs) 'ground-potential)) 

or with the use of a ground predicate: 

(defpred ser ies-1-topology (cs rl r2) 
(node Mtl ,cs) Mtl ,rl)) 
(node Mt2 ,rl) Mtl ,r2)) 
(ground Mt2 ,r2) Mt2 ,cs))) 

The current constraint for the ground node is redundant, and it seems 
that it is not very useful in constraint propagation, so it is omitted. Now the 
problem of determining the voltage across the current source in the series circuit 
is solved by determining the node potentials. The determined value for this 
voltage, not surprisingly, turns out to be independent of the ground potential 
(this independence is gotten via the algebraic simplifier). 

In more complex circuits the choice of a reference node (ground) can be 
important to the reasoning process. The absence of an explicit KCL constraint 
for the ground node can prevent constraint propagation in some cases, in others 
the unconstrained plunk placed on the ground potential can fail to induce 
propagation. Whether or not these failures in the constraint propagation occur is 
largely dependent on which node is chosen for ground. Thus the choice of a 
ground node can greatly effect the ease with which circuits are analyzed by the 
system. 

The extra degree of freedom in the node potentials is not present in the 
branch voltages, which are differences between node potentials. All currents in 
a circuit are functions of the branch voltages of that circuit and do not depend 
on the absolute potentials of the nodes. Thus circuit analysis is ultimately 
concerned only with branch voltages and not with the absolute potentials. 
Therefore all of the quantities which are of interest in electrical analysis do not 
depend on the unconstrained potential assigned to a reference node. Therefore 
the choice of the reference node does not effect the form of expressions for 
quantities of interest. However, as will be seen below there are cases in which 
the choice the quantity to be plunked has a great effect on resulting expressions 
for quantities of interest. 

Another example of the use of unconstrained plunks is circuits with 
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external terminals. In such circuits one is interested in the constraints the circuit 
imposes on the currents and potentials of the external terminals. To put the 
discussion in a more concrete framework consider the following circuit: 



/ o- 



x o- 





Figure 5. The "3-Paratlel" Circuit 

(defpred 3-parallel (circuit) 
(one-port circuit) 
(resistor ' (rl .circuit)) 
(resistor ' (r2 .circuit)) 
(resistor ' (r3 .circuit)) 
(term-eq '(tl .circuit) 

'(composite (» tl rl .circuit) 

(composite (» tl r2 .circuit) 

(» tl r3 .circuit)))) 
(term-eq *(t2 .circuit) 

'(composite (» t2 rl .circuit) 

(composite (» t2 r2 .circuit) 

(» t2 r3 .circuit))))) 

None of the branch voltages or branch currents can be determined if the 
environment of the circuit is left unspecified. However it is possible to reason 
about this circuit by giving the voltage across the terminals of the circuit an 
unconstrained plunk. The following dialog with the system demonstrates the use 
of the unconstrained plunk: 
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(assert (3-parallel 'cl)) 
(3-PARALLEL Cl) 

(what-is '(current cl)) 
(CURRENT Cl) 

(assert (== '(voltage cl) 'input-voltage)) 
(== (VOLTAGE Cl) INPUT-VOLTAGE) 

(uhat-is '(current cl)) 
(// (* INPUT-VOLTAGE 

(+ (* (RESISTANCE (Rl CD) 

(+ (RESISTANCE (R2 CD) 
(RESISTANCE (R3 CD))) 
(* (RESISTANCE (R2 CD) 
(RESISTANCE (R3 CD)))) 
(* (RESISTANCE (Rl CD) (RESISTANCE (R2CD) (RESISTANCE (R3C1)))) 

The constraint propagation to produce this result is straightforward. 
The substitution process ensures that the voltage across the external terminals of 
the circuit is in the same equivalence class as the voltage across each of the 
resistors. Thus when the external voltage is plunked the voltage across each 
resistor is determined and therefore the current through each resistor is 
determined via Ohm's law. The current through the external terminals is then 
determined as a function of the voltage across them. Unfortunately the system is 
not capable of generalizing such knowledge and applying it to situations in which 
this circuit is embedded in larger circuits. However a human engineer could use 
the system to analyze general instances of circuit fragments and then incorporate 
the results into TMS predicates for those fragments. 

Other quantities might have been plunked in the above circuit. For 
example the following dialogue shows how the voltage could be solved for in 
terms of the current: 

(assert (3-parallel 'c2D 
(3-PARALLEL C2) 

(assert (== '(current c2) 'input-current)) 
(« (CURRENT C2) INPUT-CURRENT) 
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(what-is * (vol tage c2)) 
(VOLTAGE C2) 

(solve-for '(voltage c2)) 
(// (* INPUT-CURRENT 

(RESISTANCE (Rl CD) 
(RESISTANCE (R2 CD) 
(RESISTANCE (R3 CD)) 
(+ (* (RESISTANCE (Rl CD) 

(+ (RESISTANCE (R2 CD) 
(RESISTANCE (R3 CD))) 
(* (RESISTANCE (R2 CD ) 
(RESISTANCE (R3 CD)))) 

Normal plunking had to be initiated via solve-for to perform the 
analysis but the desired result was obtained. The circuit can be viewed as 
providing a constraint on its terminals which allows the current to be derived 
from the voltage or the voltage from the current. However in each of the above 
analyses, one of these quantities was chosen as the "input" and the other was 
derived as an expression in terms of it. In more complex circuits particular 
terminals can be specified as input and output terminals. One is usually only 
interested in expressions for the output quantities in terms of the inputs. Thus 
the multidirectional view of a circuit as a constraint is replaced by a 
unidirectional relationship between inputs and outputs. 

Slices 

This section develops the use of "slices" in circuit analysis. The term 
slice was coined by Gerald Sussman to refer to a form of multiple description 
which is useful in constraint propagation [Sussman 77] [Steele & Sussman 79]. 
The basic technique is to state an equivalence between a part of the structure 
being analyzed and a different structure which exhibits the same behavior. Such 
equivalences are quite naturally stated in the equality system. The simplest use 
of slices in electrical analysis is in series-parallel reduction. Consider two resistors 
in series connected to a voltage source as is shown in figure 6. 
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Figure 6. The "Series-2" Circuit. 

(defpred series-2 (circuit) 

(voltage-source ' (vs .circuit)) 

(resistor ' (rl .circuit)) 

(resistor ' (r2 .circuit)) 

(node '(>> tl vs .circuit) ' (» tl rl .circuit)) 

(node '(>> t2 rl .circuit) ' (» tl r2 .circuit)) 

(ground ' (» t2 r2 .circuit) ' (» t2 vs .circuit))) 

Now if it is asserted that some circuit, say c2, is a series-2, and then the 
system is asked for the current of rl, the system will simply respond with (» 
current rl c2). In other words the circuit will not be analyzed. The potential of 
the top node can be determined to be the ground potential plus the strength of 
the source. However the potential of the node connecting rl and r2 does not 
take on a determined value since the neither resistor has a determined voltage. 
There is no way to derive the resistor voltages since the currents are not known. 

This circuit can be solved however with the use of slices. An slice is an 
alternate description of some portion of a structured object. In circuit analysis 
slices are installed by giving designators for terminal currents and potentials in 
terms of equivalent circuits. In this case rl and r2 are equivalent to a single 
resistor with a resistance equal to the sum of their resistances. Figure 7 shows 
this circuit with the slice imposed. 
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Figure 7. A Slice on the "Series-2" Circuit. 

This equivalence is stated by defining a resistor with the appropriate 
resistance and then stating equalities between its terminal potentials and currents, 
and terminal potentials and currents in the circuit. This can be done in the 
equality system as follows: 

(defpred term-eq (tl t2) 

(== '(potential ,tl) '(potential , t2) ) 
(== '(current ,tl) '(current , t2))) 

(assert (series-2 'cl)) 

(assert (c« '(resistance series-r) 

' (+ (» resistance rl cl) (» resistance r2 cl)))) 

(assert (term-eq ' (» tl rl cl) ' (tl series-r))) 



(assert (term-eq ' (» t2 r2 cl) ' (t2 series-r))) 

With such constraints installed the current in the equivalent resistor is 
determined from the strength of the voltage source and the resistance of that 
resistor. This current then determines the current in rl and r2, which determines 
the voltages across them and therefore the potential of the central node. 

The following predicate and operator definitions give a convenient means 
of creating series and parallel slices. 

(defpred one-port-eq (tl t2 one-port) 
(term-eq tl '(tl , one-port)) 
(term-eq t2 '(t2 ,one-port))) 
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(evalfun series-eq (rl r2) 

(-> (predand (resistor rl) (resistor r2)) 

(predand (resistor '(series-resistor ,rl ,r2)) 

(c=- '(resistance (series-resistor ,rl,r2)) 

'(+ (resistance ,rl) (resistance ,r2)))))) 

(evalfun parallel-eq (rl r2) 

(-> (predand (resistor rl) (resistor r2)) 

(predand (resistor '(parallel-resistor ,rl ,r2)) 

(c== '(// 1 (resistance (parallel-resistor ,rl,r2))) 
'(+ (// 1 (resistance ,rl)) 

(// 1 (resistance ,r2))))))) 



Of course more complex circuits can be defined, and the machinery has 
now been developed to do series-parallel reduction. Consider the ladder network 
of resistors shown in figure 8. 
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Figure 8. The "Ladder" Circuit. 
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(defpred ladder (c) 

(voltage-source ' (vs ,c)) 

(resistor ' (rl ,c)) 

(resistor ' (r2 ,c)) 

(resistor ' (r3 ,c)) 

(resistor * (r4 ,c) ) 

(ladder-topology ' (vs ,c) '(rl ,c) * (r2 ,c) * (r3 ,c) ' (r4 ,c))) 

(defpred ladder-topology (vs rl r2 r3 r4) 
(node *(tl ,vs) '(tl ,rl)) 
(node *(t2 ,rl) '(tl ,r2) *(tl ,r3)) 
(node Mt2 ,r3) Mtl ,r4)) 
(ground Mt2 ,r4) Mt2 ,r2) Mt2 ,vs)) 
( ladder-topology2 rl r2 r3 r4)) 

(defpred ladder-topology2 (rl r2 r3 r4) 

(one-port-eq '(tl ,r3) Mt2 ,r4) Mseries-eq ,r3 ,r4)) 
(one-port-eq '(composite (tl ,r2) (tl (series-eq ,r3 ,r4))) 
'(composite (t2 ,r2) (t2 (series-eq ,r3 ,r4))) 
' (paral lel-eq ,r2 (series-eq ,r3 ,r4))) 
(one-port-eq '(tl ,rl) 

*(t2 (paral lel-eq ,r2 (series-eq ,r3 ,r4))) 

'(series-eq ,rl (paral lel-eq ,r2 (series-eq ,r3 ,r4))))) 

Here is a dialog with the system about this circuit: 

(assert (ladder *c2)) 

T 

(uhat-is ' (>> current r4 c2)) 

(// (* (» STRENGTH VS C2) 

(RESISTANCE (PARALLEL-EQ (R2 C2) (SERIES-EQ (R3 C2) (R4 C2) ) ) ) ) 
(* (RESISTANCE (SERIES-EQ (Rl C2) 

(PARALLEL-EQ (R2 C2) 

(SERIES-EQ (R3 C2) (R4C2))))) 
(RESISTANCE (SERIES-EQ (R3 C2) (R4C2))))) 

This value was derived by straightforward propagation. The potential 
difference across the equivalent resistor for the entire circuit is determined by the 
voltage source. The current through that resistor is then determined by Ohm's 
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law. This determines the current through rl and through the equivalent 
resistance for the remainder of the ladder. The current through this equivalent 
resistor determines the potential across it via Ohm's law, which in turn determine 
the potential across the pair of resistors r3 and r4. The current through this pair 
of resistors is determined via the use of their equivalent resistance and Ohm's 
law. This finally determines the current in r4. 

The value given for the current in r4 is in terms of the equivalent 
resistances used in the slices. It is not immediately clear whether the better-name 
relation should consider this a better or worse representation than an algebraic 
expression in terms of the resistances of the actual resistors in the circuit. It can 
be argued that the use of the equivalent resistances prevents the algebraic 
simplification of expressions. However in the above case the expression 
containing the equivalent resistances is actually shorter (by 6 symbols) than the 
result of simplifying an algebraic expression in the original resistances, which is 
shown below: 

(// (* (+ (» RESISTANCE R3 C2) (» RESISTANCE R4 C2)) 
(» STRENGTH VS C2)) 
(+ (* (» RESISTANCE Rl C2) 

(+ (» RESISTANCE R2 C2) 
(+ (» RESISTANCE R3 C2) 
(» RESISTANCE R4 C2)))) 
(* (» RESISTANCE R2 C2) 
(+ (» RESISTANCE R3 C2) 
(» RESISTANCE R4 C2))))) 

Not enough experience has yet been had with the electronics system to 
tell whether the use of the equivalences in the results of analysis is always a good 
thing. The present system uses them. 
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"Redrawing" Circuits 

The equivalences used to solve the above ladder circuit do not always 
result in a complete circuit analysis. There are cases in which series parallel 
reduction should clearly be possible, but that the slices given so far are not 
adequate. Consider the circuit shown in figure 9. 




Figure 9. Two Sets of Parallel Resistors in Series. 

A parallel series reduction of this circuit can be attempted by placing 
parallel slices on the two pairs of parallel resistors in the circuit. Since the 
strength of the current source is a transparent designator, propagation begins 
from that point. The sum of the currents into the upper terminals of rl and r2 
is determined via the current constraint at the top node. This determines the 
current in the parallel equivalent of those resistors, which determines the sum of 
the currents in the lower terminals of rl and r2. Now it would seem that the 
sum of the current in the upper terminals of r3 and r4 should be determined. 
However if the node predicate has been used in the standard way to create the 
current constraint for the center node, then the constraint net does not constrain 
the sum of the currents in r3 and r4 in a manner which allows this determination 
to be used. The current constraint for the node connecting rl, r2, r3 and r4 
might look like: 
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(== ' (>> current t2 rl) 

* (- (+ (>> current t2 r2) (» current tl r3) (» current tl r4)))) 
(== ' (>> current t2 r2) 

* (- (+ (» current t2 rl) {» current tl r3) (» current tl r4)))) 
(== ' (>> current tl r3) 

' (- (+ (» current t2 rl) (» current t2 r2) (» current tl r4)))) 
(== ' (>> current tl r4) 

' (- (+ (» current t2 rl) (» current t2 r2) (» current tl r3)))) 

This constraint cannot give a better name to the sum of the currents in 
the upper terminals of r3 and r4 simply because none of the equalities 
representing the constraint deal with that sum. This prevents the current in the 
equivalent resistance of r3 and r4 from being determined. There is no other way 
for constraints to propagate around the circuit. Even though the voltage across 
rl and r2 can be determined, and therefore the two branch currents, the 
potential difference across the current source has not been determined, so that 
the potential across the bottom resistors remains unknown. Also, since no 
current constraint is created at the ground node, the sum of the currents in r3 
and r4 are not directly constrained. (If such a current constraint were present an 
example of the failure of parallel slices would only be slightly more complex.) A 
solution to this dilemma is suggested by the alternate drawing of the circuit 
shown in figure 10. 

The parallel slices are somehow more strongly suggested in this drawing. 
This is due to the explicit terminals which can then be made equivalent to the 
terminals of the resistors in the slices. As a result of this, the current constraints 
implicit in the diagram more directly constrain the current in the equivalent 
resistance. It is possible to "redraw" a circuit by restructuring the current 
constraints to take into account terminals of slices. This can be done directly 
using the predicates and operators defined so far. For example a parallel-series 
reduction of the above circuit can be accomplished if the following is used to 
construct the central node: 

(connected '(composite (t2 rl) ' (t2 r2)) '(composite (tl r3) (tl r4))) 

The present method of creating slices relies heavily on the user of the 
system. It would be far more satisfactory to have the system recognize the 
appropriate equivalences on its own. This has not been accomplished to date and 
might represent a fruitful area for further research. 



Chapter 



44 



Electronics 




Figure 10. An Alternative Drawing. 

The major point brought out in the discussion of electronics presented in 
this chapter is that an effective reasoning system can be derived from algorithms 
designed to handle the problems surrounding equality and multiple descriptions. 
The majority of the reasoning strategies used in these algorithms can be justified 
in terms of general principles without reference to electronics. 
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Chapter IV 
Algorithmic Details and Relation to Other Work 

The Truth Maintenance System 

A truth maintenance system or TMS is used in the equality system to 
record and enforce logical relations among propositions. The basic functions 
performed by the TMS were first introduced by Richard Stallman and Gerald 
Sussman in an electrical analysis system which made assumptions about transistor 
states [Stallman & Sussman 77]. A separate module, called a truth maintenance 
system, or TMS, was later designed by Jon Doyle [Doyle 78]. Doyle's TMS keeps 
track of propositional justifications and uses them to incrementally update beliefs 
and track down assumptions underlying contradictions. The propositional 
reasoner described here is a further refinement of these ideas and is very similar 
to a truth maintenance system developed by the author and described in a 
separate publication [McAllester 78]. The basic principles of that TMS are 
briefly stated here and some comments about its relation to other are given. 

Each assertion or belief in the system is given a TMS node which can 
take on one of three truth values, true, false, or unknown. Logical relations 
among beliefs are recorded in disjunctive clauses which should be viewed as 
simple constraints on the truth values associated with nodes. For example the 
representation of an implication between a tms node nl, and a second tms node 
n2 would be (or (nl . false) (n2 . true)). If the assertions represented by nl 
and n2 were mutually contradictory (could not both be true at the same time), 
this would be represented in the TMS by the clause (or (nl . false) (n2 . 
false)). In general a clause is a list of terms, each of which is an association of a 
TMS node with either "true" or "false". A clause states a constraint on the truth 
values of the TMS nodes which says that one of the nodes in the clause must 
have the associated truth value. 

A clause can locally determine the truth value of a node. Whenever all 
the nodes in a clause except one have the opposite truth value from the one 
associated with it in the clause, the clause has only one chance left to be 
satisfied. When this occurs, and the truth value of the node which might still 
satisfy the clause is "unknown", the truth value of that node is set to the value 
associated with it in the clause. A pointer is constructed from the node whose 
truth value was deduced to the clause which deduced it. This pointer gives the 
well founded support for the truth value of the node and is used in generating 
explanations and during dependency directed backtracking. 

The user of the TMS can add TMS clauses at will, but once in the 
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system they cannot be removed. This poses absolutely no problem since 
removable clauses can always be simulated by creating a TMS node to represent 
the truth of the clause. This is done by the -> predicate which creates a TMS 
node to represent an implication and adds the clause stating that the implication 
along with the implication's antecedents implies the conclusion. The TMS node 
representing the implication is then returned to the user who can set its truth 
value as he pleases, effectively adding and removing the implication represented 
by that node. 

When the user sets the truth of a node to "true" or "false", as can be 
done with assert or set-truth, this truth value is taken by the system to be a 
premise. As premises are added the TMS enforces the constraints represented by 
its clauses and deduces other truth values. The user can also retract premises via 
a remove-truth function which makes a node unknown. When this is done the 
support for other truth values can become invalid. Whenever the supporting 
clause for the truth value of a node can no longer be used to deduce the value of 
that node the supported truth value is removed. Thus the removal of a truth 
value can propagate to generate the removal of a large number of other truth 
values. After this propagation has occurred each of the nodes whose truth value 
was removed must be checked for alternate supports for truth values. By doing 
this only after the removal phase is completed the system avoids the possibility of 
looping support structures. 

Because constraints can be added at any time, and because loops can 
exist in the constraint set, it is possible for contradictions to arise. A 
contradiction is nothing more than a clause in which all the nodes have the 
opposite truth value from the one associated with it in the clause. When this 
happens the support pointers of the nodes in the clause can be used to track 
down the set of premises which underlie the truth values of those nodes. One of 
these premises can then be retracted to remove the contradiction. 

Whenever a truth value for a node is an underlying support for the 
truth value of a node in a contradiction, i.e. whenever a truth value leads to a 
contradiction, it is valid to deduce the opposite truth value for that node. There 
may be loops in the constraints in the TMS which prevent this from happening 
naturally. For this reason a backtracking facility is provided which will take a 
contradiction and an underlying premise and use clause resolution to generate 
new clauses which bypass the loops preventing the deduction of the negation of 
the premise. Since small clauses are useful in more situations, the backtracker 
does not compute a "conditional proof" constraint involving all the premises 
underlying a contradiction as is done in other systems [Doyle 78]. It instead does 
a minimum amount of clause resolution which adds clauses just large enough to 



Chapter IV 47 Algorithmic Details 



break the constraint loops. 

Demons can be attached to the TMS nodes which are invoked when the 
node changes truth value. There are three types of demons which can be 
attached to a node, true-noticers, false-noticers, and unknown-noticers, which are 
activated when the node takes on the suggested value. True-noticer demons are 
attached to equality nodes in the equality system which activate certain canonical 
naming functions when the equality becomes true. Nodes can be given unknown- 
noticer demons which set the node to a default truth value (as a premise) when 
the node would otherwise be unknown. These noticers are run only at times 
when the TMS is otherwise stable, making interactions between TMS functions 
and noticer actions comprehensible. 

Assumptions are simply specially marked premises. This marking serves 
no other purpose than to act as a guide when choosing a premise to retract 
during backtracking. There is no significant distinction in this system 
between an assumption and a simple premise. Non-monotonic dependency 
structures have been used in other systems to make assumptions by justifying a 
belief in terms of a lack of knowledge to the contrary [Doyle 78]. Such 
dependency structures were found to be totally unnecessary in the current system 
and in the author's opinion they lead to considerable algorithmic complexity and 
conceptual obscurity. With a slight modification to the TMS as described here it 
can also be shown that there is no expressive power gained from non-monotonic 
dependencies. 

The Canonical Naming Algorithm 

The basic mechanism which must be implemented is the determination 
of a canonical name for an equivalence class under the equalities provided to the 
system. The equalities can be though of as arcs between nodes representing the 
designators in a graph. Finding a canonical name involves finding a designator 
which is at least as good as any other designator in the class under the better- 
name relation. A marking algorithm is used to scan the designators in a class. 
The marker placed on a class of designators is maintained on each designator in 
that class and a pointer is maintained from the marker to the canonical name of 
that class. This makes access to the canonical name of any designator's 
equivalence class very efficient. 

A major design goal is the ability to add and remove equalities at will, 
and maintain proper canonical names. This algorithm must therefore be extended 
to handle the maintenance of the canical names when the set underlying 
equalities changes. The freedom to add and remove equalities means that 
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equivalence classes can merge and split as the equalities change. I will first 
consider the case of adding equalities. When a new equality is added (its TMS 
node takes the value "true"), it can merge two previously distinct equivalence 
classes. Thus each new equality added must be checked to see if a marker can 
be propagated across it. If the two designators which become linked have 
distinct markers on them, then one of the markers, call it ml, is chosen to 
propagate across the equivalence class. The other marker, call it m2, can be 
removed from all designators on which it appears since ml must propagate across 
the entire class previously marked by m2. As a marker propagates across new 
members of its equivalence class the canonical name pointer must be 
appropriately kept up to date. 

Now consider the problems involved in retracting equalities. Each 
marker starts from a single designator and propagates to other names via 
equalities. The class marked by the marker is defined to be the class containing 
the marker of origin. This class can change as equalities are added and removed. 
In light of this definition, the presence of a marker on a designator can be 
considered valid only if the designator is in the same class as the designator from 
which the marker originated. Thus the marking of a designator with a marker is 
associated with a TMS node representing the equality between the marked 
designator and the origin of the marker. 

As a marker is propagated across an equality from designator 1 to 
designator a logical relation is added to the TMS which states that the equality 
between designator 1 and the marker origin, along with the equality between 
designatorl and designator, imply the equality between designator and the 
marker origin. The truth of the equality between designator and the origin is 
always checked when the marker is used to find a canonical name for 
designator. If any equality is later removed, then the equalities which depend 
upon it are removed (their TMS nodes take on a truth value of "unknown") and 
marking associated with these equalities can be recognized as invalid. 

Because equalities can be retracted at any time, it is possible that the 
canonical name pointed to by a marker could be removed from the class marked 
by the marker. Since the primary purpose of the markers on designators is to 
point to the canonical name for that designator, this would render the marker 
useless. For this reason a marker is only considered valid while the canonical 
name it points to is in its equivalence class. A canonical name function has been 
written which takes any designator, first ensures that that designator is validly 
marked by a valid marker and then returns the canonical name pointed to by 
that marker. In this way the canonical names are correctly maintained under the 
retraction of equalities. 
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Substitution 

The means of controlling the substitution process have not been 
discussed so far. An image function is defined on designators such that the 
image of a designator is the result of replacing all the top level subdesignators by 
their canonical names. The image of a name is equal to the name as long the 
equalities between the parts hold. The invariant which the equality system 
attempts to maintain is that the images of all designators are explicitly 
represented in that designators class. This allows the image of the designator to 
be a candidate for the canonical name of that class. Since canonical names are 
continuously chang..ig, the image of a designator changes also. Thus a designator 
must be monitored in some way to ensure that its image is always in its class. 

Before discussing the details of the way in which this process is 
controlled, first consider some general properties of the image function. An 
internally-equal relation can be defined on designators which is stricter than the 
standard equality relation used throughout this document. Two designators are 
internally-equal if they have the same number of top level subdesignators and 
each pair of corresponding top level subdesignators are equal. All internally-equal 
designators have the same image. For example (+ 1 2) is internally equal to (+ 
(- 2 1) (+ 1 1)) but not internally-equal to 3 or (+ 2 1), assuming the 
appropriate numerical equalities have been added. Both (+ 1 2) and (+ (- 2 1) 
(+1 1)) have the same image, which should be (+ 1 2). An equivalence class 
can be divided into a set of internally-equal subclasses. Only one element from 
each of these subclasses need be monitored with respect to substitution since all 
the elements in a given class have the same image. 

All designators have their image computed at least once. A TMS node 
is constructed to represent the assertion that the image is internally-canonical, i.e. 
all its top level subdesignators are the canonical names of their equivalence 
classes. Each canonical name has a TMS node associated with it representing the 
assertion that it is the canonical name for its equivalence class, and the truth of 
this node is maintained by the equality system. These nodes can be used to 
construct a support for the truth of the node representing the assertion that an 
image is internally canonical. A demon is placed on this internally canonical 
node which recomputes the image of that designator when the node becomes 
unknown, thus monitoring that designator in a way that ensures its image will be 
recomputed when needed. This single monitoring is all that is needed to monitor 
the entire subclass of designators which are internally-equal to this designator. 
Thus when an image is computed only the image is monitored for future need to 
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recompute an image. 

Because equalities can be retracted, it is possible that internally-equal 
subclasses get split up. When this happens the single designator which was being 
n nitored is no longer sufficient to monitor all the designators which were in 
that subclass. To make sure that all designators are properly monitored, an 
internally-equal node is associated with every designator which is being monitored 
via its internal equality to another designator. This internally-equal node 
represents the assertion that the designator is internal equal to some other 
designator which is monitored. When an image is computed from a designator, 
and is not equal to that designator, a node representing the internal-equality 
between the two designators is created and associated with the designator whose 
image was computed. If this internally-equal node ever becomes unknown, then 
the image of that designator is recomputed. Thus it is possible to ensure that the 
image of every designator is always in that designators equivalence class and has 
a chance at becoming the canonical name of that class. 

The present equality system runs unreasonibly slow, on the order of a 
minute per constraint in constraint propogation analysis of electronic circuits. 
The primary bottleneck in the process turns out to be inefficient hashing of 
designator expressions. Improving the hashing methods should yeild about an 
order of magnitude improvement in computation time, but the resulting six 
seconds per constraint is still too slow to yeild a practicle system for serious 
problems. This is a first pass implementation and it is hoped that future 
refinements can still further significantly reduce the computation times involved. 

Relation to Other Mechanisms for Handling Equality 

This section discusses several classes of work related to the equality 
system described in this document. The first relates the work discussed here to 
previous work which used canonical expressions to represent equivalence classes. 
The second discusses an approach taken to equality in the resolution theorem 
proving tradition. A third approach to equality incorporates it into pattern 
matching (or unification) procedures. Finally, special attention is paid to an 
algorithm developed by Knuth and Bendix to decide equivalences in certain 
universal algebras. 

Canonical forms have been used to represent classes of equivalent 
expressions in many previous systems. Ira Goldstein used canonical forms for 
geometric objects in a geometry theorem prover [Goldstein 73]. More recently 
Howard Shrobe used canonical representatives of equivalence classes in a more 
general purpose reasoning framework [Shrobe 79]. In Shrobe's system substitution 
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was performed in such a way that all reasoning took place in terms of the 
canonical expressions. There .are certainly many other systems which employ 
similar methods that the author is unaware of, and forgiveness must be asked 
from those who are not referenced. The originality claimed here is not the 
concept of using canonical representatives of equivalence classes but a particular 
application of this mechanism to reasoning. In the systems mentioned above the 
existence of equivalent expressions is dealt with as a problem to be overcome in 
the reasoning process. Canonical expressions are used to unify the knowledge 
about an individual expressed in terms of different designators for that expression. 
If equality is only viewed from this perspective then one tends to minimize the 
role of multiple designators for the same referent in the reasoning process. In the 
applications of equality discussed here, aspects of the canonical naming 
mechanism which are ignored by these other systems, such as substitution into 
non-canonical expressions, and the detailed structure of a better-name relation, 
take on a great significance in the reasoning process. Here equality is viewed as 
playing a central role in reasoning, rather than as a necessary evil. 

A method for handling equality, called paramodulation, has been 
developed in the context of resolution theorem proving [Robinson 65] [Chang & 
Lee 73]. This technique is a clause resolution rule which incorporates 
substitution of equals for equals. This allows deduction of anything deducible in 
a system of first order predicate calculus with equality. The paramodulation rule 
itself however provides no guidance as to when substitutions should be done and 
systems using little or no heuristic guidance become lost in combinatorial 
explosions. Several methods have been developed which cut down on the number 
of allowed applications of paramodulation, such has hyperparamodulation and 
linear paramodulation. However such restrictions on paramodulation do not give 
the detailed guidance of substitution provided by canonical naming schemes. 

Another approach taken to equality involves incorporating equality 
axioms into pattern matching (or unification) [Fay 78]. A discussion of this is 
best done by an example. Suppose that op is an associative operator (i.e. (op (op 
x y) z) = (op x (op y z))). The match of (op (op a b) c) with (op x y), 
where x and y are variables and a, b, and c are constants, would yield two 
substitutions: ((x . (op a bll , (y . c)} and {(x . a) , (y . (op b c))j. 
Thus the equality variants of the two expressions being unified are taken into 
account in generating the variable bindings. As in the case of paramodulation, 
such a scheme can be incorporated into a system that can make all valid 
deductions from a given set of axioms in first order predicate calculus. However, 
for this approach to be feasible, the equality axioms must be determined at the 
time the matching is done. This requirement makes this approach inappropriate 
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to the type of reasoning about equality presented in this document where the set 
of useful equalities is continually expanding. 

The most notable of the algorithms for handling equality is one 
developed by Donald Knuth and Peter Bendix [Knuth & Bendix 69]. They give 
an algorithm for deciding whether two words (designators) are equivalent under 
equivalence relations defined by simple sets of algebraic axioms. A partial order 
on words, analogous to the better-name relation used here, is used to reduce a 
given word to a canonical form. Given two words it is possible to tell if they are 
equivalent by seeing if their canonical forms are the same. Knuth and Bendix 
represent the algebraic axioms as a set of reductions which take the form of 
rewrite rules and represent equalities. They also give a decision procedure for 
determining if a given set of reductions is complete, i.e. results in a complete 
decision algorithm for tests of equality under those axioms. Perhaps the most 
interesting aspect of their work is a method given for extending a reduction set 
to a complete set by adding derived reductions. 

When comparing the algorithms developed by Knuth and Bendix to the 
equality system described here, a strong similarity emerges. Each equality in the 
equality system acts as a reduction in the sense that it can result in one 
designator being replaced by another inside a larger designator when substitution 
occurs. Reductions which contain variables can not be encoded as equalities in 
the equality system. However arbitrary reductions can be placed in the 
evaluation functions. Since the general equality system allowed for an arbitrary 
better name relation, the relations used by Knuth and Bendix could be easily 
incorporated. Therefore any partial order on designators and any complete 
reduction set used by a Knuth and Bendix type algorithm can be used in the 
equality system. Furthermore since each derived equality acts as a new 
reduction, any reduction set which can be extended to a complete set by the 
Knuth and Bendix algorithm will, in much the same way, automatically form a 
complete set in the equality system. 

I do not mean to claim a rediscovery of the results on completeness 
found by Knuth and Bendix which involve the developement of a specific better- 
name relation and its interaction with specific sets of equality axioms. I am not 
particularly concerned with completeness however, and approached the problems 
presented by equality from a completely different perspective. They were 
concerned with deciding equality under a small set of algebraic axioms which 
remain fixed over time. I am concerned with the use of equality as a means of 
knowledge representation and as a mechanism for handling multiple descriptions. 
Thus I am concerned with a very large rapidly expanding set of equalities (or 
reductions) which can be manipulated as the reasoning process progresses. 
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An important fundamental distinction between the system described here 
and all previous work on equality (at least to the author's knowledge) is that this 
system views equality and multiple descriptions as a central part of the reasoning 
system. Thus it makes sense to ask the system for the value of an expression 
even though the expression itself is a designator for that value. In traditional 
theorem proving systems such a request could not even be properly stated, an 
expression for the desired quantity is trivially found and the system has no way 
of telling that this does not completely solve the problem at hand. 
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Appendix One 

A Basic TMS and Equality System User's Manual 

The basic utility functions of the TMS and the equality system are 
divided into three sections. First are the TMS utilities which are used in the 
code for the equality system. These primitives allow the direct creation of TMS 
clauses and the manipulation of the truth values associated with TMS nodes. 
However they are not intended to be used at the interactive user level. The next 
set of functions are TMS utilities which are intended to be used interactively by 
a system user, including the basic TMS predicates. The final section presents the 
basic equality system utilities. 

Internal TMS Utilities 
DEPENDENCY-NODE 

(dependency-node osser t ion>) 

This function takes an "assertion" in the form of an s-exprcssion 
and returns a TMS node which has been associated with that assertion. 
The assertions are kept in a hash table which is searched each time 
dependency-node is called. Thus repeated calls with the same assertion 
yield the same node. 

SET-TRUTH 

(set-truth <node> <value> <setter>) 

This function takes a node, a value which is either the atom true 
or the atom false, and an atom used to record the source of the set value. 
The node takes on the truth value and- the TMS automatically propagates 
the results of that value. Noticer functions, which are meant to notice 
changes in truth values, are called automatically when the propagation has 
stabilized. A noticer function can induce further changes in truth values 
and the process continues until no noticers are left to run. 

The most common value for the setter argument is the atom 
"premise". However, if the setter argument is the atom "default", then the 
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set truth value is marked as an assumption and is treated specially during 
backtracking. 

If the node already has a truth value, the node is still given the 
new value with the new justification. This can lead to contradictions if 
the opposite value had been previously deduced. Such contradictions are 
treated the same as any other contradictions in the system. 

REMOVE-TRUTH 

(remove-truth <node>) 

This function is used to retract premises. The truth value of 
<node> will be retracted and the appropriate propagation of truth value 
removal and invocation of noticer functions is done. 

T-NOT 

(t-not <term>) 

A term is either a node or an association of a node with either 
the atom true or the atom false. It is often convenient to be able to 
bundle an node with a truth value in using primitive TMS functions. The 
t-not primitive returns an association of a node with the opposite value 
from the one associated with it in the argument to t-not. A single node is 
always treated as if it were associated with "true". 

IMPLIES 

(implies <terms> <term>) 

This function takes a list of antecedent terms and a consequent 
term and installs the constraint on node truth values corresponding to this 
implication. Once the constraint has been added it cannot be removed. 
For this reason the primitive implies should only be used to state 
tautologies, or relations which are true by definition. 
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CONTRADICTORY 

(contradictory <terms>) 

This primitive is essentially the same as implies. It takes a list of 
terms and installs the constraint that at least one of the nodes in those 
terms must have the opposite value from the one associated with it in the 
term. Again this should be used only to state tautologies. 

TRUE-NOTICER, FALSE-NOTICER, UNKNOWN-NOTICER 

( true-not i cer <node> <form>) 

These functions attach noticers to nodes. The noticer is a form 
to be evaluated when the node takes on the appropriate truth value. 

DEFAULT 

(default <node> <value>) 

If node has an unknown truth value then this sets the value to 
the one given with a setter of "default", which marks the value as an 
assumption to be treated specially by the backtracker. It also places an 
unknown noticer on the node which will set the node to the default value 
whenever it would otherwise be unknown. 

IMPLIES-UM 

(implies-uni <terms> <term>) 

This is just like implies except that the constraint generated will 
not automatically deduce the negation of any antecedent term. If the 
constraint is violated however, it is treated as a contradiction and can then 
be used to deduce the negation of any antecedent term. This was 
implemented primarily to capture the full expressive power found in other 
systems which use non-monotonic dependency structures. The actual 
utility of this feature is an open question. 
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User Level TMS Utilities 

These utilities are intended to be used at the interactive user level. To 
achieve sufficient generality they take "propositions" as arguments. A 
"proposition" is either a term (a TMS node or a TMS node associated with a 
truth value) or an s-expression representing an assertion which is coerced to a 
node via dependency-node as described in the previous section. 

ASSERT 

(assert <prop>) 

This function takes a proposition and calls set-truth on the 
corresponding node. The atom premise is always used as the setter. Thus 
(assert <node>) is identical to (set-truth <node> 'true 'premise). And 
(assert 'p) is identical to (set-truth (dependency-node 'p) 'true 

'premise). 

RETRACT 

(retract <prop>) 

If the node associated with <prop> has a truth value which is 
given as a premise then this function removes that truth value values (and 
all values which critically depended on that value). If the proposition has 
a deduced truth value or a default truth value then this function is 
effectively a no-op. 

ASSUME 

(assume <prop>) 

This sets the truth value of the node associated with the 
proposition to a truth value and marks that node as an assumption which 
can be recognized by the backtracker. It does this via a call to defaul t as 
described above. 
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UNLESS 

(unless <propl> <prop2>) 

This gives prop2 a true default value and installs a unidirectional 
implication such that if propl becomes true prop2 will be deduced to be 
false. However since this implication is unidirectional the assumption that 
propz is true will not result in a deduction that propl is false. This 
primitive depends on the impiies-uni facility described above and is of 
questionable utility. 

PREDNOT 

(prednot <prop>) 

If the proposition is a TMS node or an assertion s-expression this 
returns a dotted pair of the TMS node and "false". Thus (prednot 
<node>) is equivalent to (cons <node> 'false). And (assert (prednot 
<node>)) is equivalent to (set-truth <node> 'false 'premise). If the 
argument to prednot is already an association of a node and a value then a 
pair of the node and the opposite value is returned. 

WHY 

(why <query>) 

A query is either a proposition or a numerical reference to nodes 
mentioned in previous explanations generated via why. Applications of this 
function can be divided into three cases. First, if the node referenced in 
the query has an unknown truth value then a simple statement to that 
effect is printed. Second, if the node has a truth value which was given as 
a premise then the truth value is stated and the user is told that it is a 
premise and gives him the value of the setter argument used when the 
premise was made. Finally, if the node has a deduced truth value then the 
associated truth value is stated along with a numbered list of supporting 
terms. The number associated with each term can be used to reference 
that term in the next use of the why function. 
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-> 



The number can be used to "pop" back to the. previous 
justification and allow numerical references to terms in that justification. 
Therefore it is very easy to walk across the support tree for any assertion 
via numerical references to the supporting terms. 



(-> <propl> <prop2>) 

If the propositions are simple nodes or assertions the above 
expression returns a TMS node which has the following TMS constraint 
imposed: 

(or ( (-> <nodel> <node2>) . false) 
(<nodel> . false) 
(<node2> . true)) 

If the propositions are simple nodes or assertions an appropriate constraint 
is installed for the associated truth values. 

PREDOR 

(predor <propl> <prop2> <prop3> ...) 

If the arguments are simple nodes or assertions the above form 
returns a node with the below TMS constraints imposed. Otherwise 
analogous constraints are imposed for the associated truth values. 



Primitive Functions 60 Primitive Functions 



(or ( (predor <nodel> <node2> ...) . false) 

(<nodel> . true) 

<<node2> . true) 
...) 

(or (<nodel> . false) ((predor <nodel> <node2> ...) . true)) 

(or (<node2> . false) ((predor <nodel> <node2> ...) . true)) 

PREDAND 

(predand <propl> <prop2> <prop3> ...) 

This is analogous to predor. 

PRESUMABLY 

(presumably <prop>) 

This predicate performs the following side effect involving the 
node returned: 

(assume (-> '(presumably <prop>) <prop>)) 
DEFPRED 

The defpred construct is a means of defining TMS predicates. A 
defpred can have the same syntax as a LISP defun. The node returned by 
a predicate defined with defpred represents the assertion that the predicate 
is true of the objects to which it was applied. The body of a defpred is a 
list of forms which evaluate to either a TMS node or nil. During an 
application of the defined predicate, if a form in the body evaluates to a 
TMS node then that node will be implied by the node returned as the 
value of the predicate application. If the form returns nil it has no effect 
on the node returned from the predicate application. Some examples ^re 
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given below: 

(defpred vertebrate (x)) 
(defpred flys (x) ) 

(defpred bird (x) 

(vertebrate x) 
(flys x)) 

(defpred hairy (x) ) 

(defpred female (x)) 

(defpred bears- I i ve-young (x)) 

(defpred mammal (x) 

(vertebrate x) 
(presumably (hairy x)) 
(presumably (-> (female x) ! 

(bears- I i ve-young x)))) 

(defpred duck-bi I I (x) 
(mammal x) 
(prednot (bears- 1 i ve-young x) ) ) 

The following is a dialogue with the system using the above predicates. 

(assert (mammal 'joe)) 
(MAflflAL JOE) 

(assert (female ' joe) ) 
(FEMALE JOE) 

(why (bears- I i ve-young 'joe)) 

((BEARS-LI VE-YOUNG JOE) IS TRUE FROM) 

(1 (FEMALE JOE) IS TRUE) 

(2 (-> (FEMALE JOE) (BEARS-LIVE-YOUNG JOE)) IS TRUE) 

T . 
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(why 2) 

((-> (FEMALE JOE) (BEARS-LIVE-YOUNG J0E>) IS TRUE FROM) 
(1 (PRESUMABLY (-> (FEMALE JOE) (BEARS-LIVE- YOUNG JOE))) IS TRUE) 
(2 (-> (PRESUMABLY (-> (FEMALE JOE) (BEARS-LIVE-YOUNG JOE))) 
(-> (FEMALE JOE) (BEARS-LIVE-YOUNG JOE))) 
IS TRUE) 
T 

(why 2) 

((PRESUMABLY (-> (FEMALE JOE) (BEARS -LIVE -YOUNG JOE))) IS TRUE FROM) 

(1 (MAMMAL JOE) IS TRUE) 

(2 (-> (MAMMAL JOE) 

(PRESUMABLY (-> (FEMALE JOE) (BEARS-LIVE- YOUNG JOE) ))) ) 
T 

(uhy 1) 

((MAMMAL JOE) IS TRUE AS A PREMISE SET BY PREMISE) 

(assert (cluck-bi I I ' joe) ) 
(DUCK-BILL JOE) 

(why (bears- I ive-young 'joe)) 
((BEARS-LIVE- YOUNG JOE) IS FALSE FROM) 
(1 (DUCK-BILL JOE) IS TRUE) 

(2 (-> (DUCK-BILL JOE) (NOT (BEARS-LIVE-YOUNG JOE))) IS TRUE) 
T 

The repeated application of a predicate to the same arguments 
returns the same node, but the body of the defpred is only evaluated on 
the first call. Thus predicate applications can be used in such places as the 
antecedents of implications, or simply in getting hold of a TMS node to 
use in a query, without of generating multiple copies of TMS constraints. 

There is an alternate syntax for defpred in which the bound 
variable list is replaced by an atom which is bound to a list of the values 
of the arguments in a application of the defined predicate. The following 
use of this syntax is from electrical circuit analysis: 
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Naming System Utilities 



'(== <designatorl> <designator2>) 

This is a predicate which returns a TMS node representing the 
equality of two designators. The node returned has a true-noticer attached 
to it that wakes up certain naming routines when the equality node 
becomes true. 

WHAT-IS 

(what-is <designator>) 

This returns the canonical name of <designator>'s equivalence 

class. 

SOLVE-FOR 

(solve-for <designator>) 

This function "solves for" a particular quantity. This is done by 
plunking quantities it can be expressed in terms of. The value found for 
the designator in question may contain plunks. If this happens the so I ve- 
for function can be applied to these plunks which initiates further 
plunking. If the constraints in the system do not determine the designator 
whose value is sought then its canonical name will be left in terms of 
plunks. 

EVALFUN 

Evalfun is used to define evaluation functions for operators which 
are used in designators. This evaluation function can be applied to any 
designator which has the appropriate top level operator. Evalfun has the 
exact same syntax as defpred. The only difference between an application 
of a predicate defined via defpred and an application of an evaluation 
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function defined by evalfua, is that terms returned by forms in the body 
of an evaluation function are made true at premises instead of being 
implied by some other node. The foflowog k aa eaampte of the use of 
evalfun for doing the addition of anaahif i with an e*a*uetioa function for 
a two place addition operator: 

(eval fun + args 

(if (apply and (mapcar 'nuaberp area)) 
(« '(4 .aargs) (apply plus args)))) 
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Appendix Two 

The Equality System Code 

This appendix contains the LISP code for the basic equality system. 
There are functions and macros used, but not defined here, which are not 
standard lisp primitives. Some of these are adopted from the LISP MACHINE 
project at the MIT AI lab. Others are parts of the authors personal utility 
package. It is hoped that the reader can infer the actions taken by these 
functions and macros from their context and pnuemonic names. However if the 
reader feels that some ambiguity exists a description of these functions and 
macros is supplied in appendix five. 
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882 (declare (special moticers* *removed-l isl* ^contradict ions* vtrace ctrace ntrace 

883 ttrue-nodef ^designators* *think-queue* *name-queue* 

884 thmk-trace 

885 . ^better-pred* *co incidence*) ) 
886 

887 (deftype (designator) 

888 desig-hash- val 

889 expansion 

818 eq-links 
811 marker 

612 eq-or iqin-node 

613 trans lat ion 
814 sub-memo-node 
615 memoirat ions 

816 st imulate-l ist 

817 uniqueness 
618 opaqueness 

819 Height 

828 plunk -weight) 

821 

822. (defmacro primitive' (name) 

823 '(atom (expansion .name))) 

624 

825 (defmacro number - ' (namo) 

826 ' (number p (expansion ,name))) 
827 

828 (defstruct (marker) 

829 origin 
838 canonical 

831 canonical -node 

832 act i ve-node) 
833 

834 (delun names-init (better-pred coincidence) 

835 (tms-init) 

836 (hash-array '^designators* 2847) 

837 (name-type 'designator) 

838 (name-type 'marker) 

839 (setq think-trace nil) 
648 (setq ttli ink -queues nil) 

841 (setq tnamn-queiie-s nil) 

842 (setq -rbet ter-predi better-pred) 

843 (setq tcoinc idence* coincidence) 

844 (setq *true-node* (dependency-node ' true-by-del ini t ion)) 

845 (set-truth strue-node* 'true 'definition)) 
846 
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882 (defun designator (exp) 

863 (designator2 (>>expand exp))) 

884 

085 (defun designator? (exp) 

886 (if (designator^ ?vp) 

887 exp 

008 (let* < (des-iqnator-exp (it (atom exp) exp (mapcar 'designator2 exp))) 

889 *as (table-asscr designator-exp '^designators*) )) 

818 (it (cdr as) 

Oil (cdr as) 

012 (let ( (n-des ignator (make-designator))) 

013 (rplacd as n-des ignator ) 

81* (sett (expansion n-designator) designator-exp) 

815 (sett (translat ion n-designator) 

816 (i I (atom exp) 

017 exp 

018 (mapcar '(lambda (des) (translat ion des) ) 

019 designator-exp))) 

020 (name- inspect n-designator) 

021 n-designator))))) 
022 

023 (delun name- inspect (name) . 

824 (II (not (primitive' name)) 

825 (let ((op (expansion (car (expansion name))))) 

026 ( i f (atom op) 

027 (mapc '(lambda (oper) 

828 (addf '(,oper '.name) 1st imulate-l ist name))) 

829 (operator-defs op)))))) 
030 

031 (defun >>expand (exp) 

032 (cond ((or (atom exp) (designator? exp)) 

033 exp) 

034 <<eq (car exp) •>>) 

035 (>>expand-2 (cdr exp))) 

036 (t (cons (>>expand (car exp)) 

037 (>>expand (cdr exp)))))> 
038 

039 (defun >>expand-2 (exp) 

040 (cond ((null (cdr exp)) 

041 (>>expand (car exp))) 

042 (t (list (>>expand (car exp)) 

843 (»expand-2 (cdr exp)))))) 

044 
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082 (defun marked' (des) 

883 (let ((ml (marker ties))) 

884 (and ml 

805 (true' (active-node ml)) 

886 (true" 1 (eq-or iqin-node des>) 

887 (let Ucan (canonical ml))) 

888 (and <<>q ml (marker can)) 

889 --(true' (eq-origin-node can))))))) 
818 

811 (detmacro insure-marked (des) 

812 '(if (not (marked' ,des>) (initiate-marker ,des))) 
813 

814 (defun class-name (des) 

815 (insure-marked dps) 

816 (canonical (marker des)) 
817 

818 (defun uhat-is <e-p> 

819 (let <(des (designator exp)>) 
828 ( insure- marked des) 

821 (think) 

022 (translation (c lass-name des) )) ) 

023 

82* (defun insure-al I -marked (desiqs) 

025 (do ((rest desiqs (cdr rest))) 

026 ((cond ((nul I rest) t) 

827 ((not (marked' (car desiqs))) 

828 (initiate-marker (car desiqs)) 

829 (insure-all-marked desiqs) 
838 t)))>> 

831 

032 (defun Image (des) 

033 (cond ((primitive' des) des) 

03* • (t (insure-all-marked (expansion des)) 

035 (let (<exp (expansion des))) 

836 (let ((im-exp (mapcar 'class-name exp) )) 

037 (if (equal im-exp exp) des (desiqnator im-exp))))))) 

038 
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662 (defun canonical' (des) 

663 ( tnsure-marked des) 

88* (let* ((ml (marker des)> 

005 (nodel (eq-or igin-node des)) 

606 (node2 (active-nod* ml)) 

007 <node3 (canonical-node ml))) 

008 (if (and U-rue' nodel) 

009 _.<tru«' node?) 

010 teq des (canonical ml))) 

011 (list nodel node2 node3)))) 
012 

013 (defun intern-canonical' (des) 
01* ( I f (primi t i ve" 1 des) 

015 (list ? true-node:) 

016 (do ((de^iqs (expansion des) (cdr desigsl) 

017 (nodes nil)) 

018 ((null desiqs) nodes) 

019 (let ((nodes: (canonical' (car desigsl))) 

020 (if nodes2 

021 (s-etq nodes (nconc nodes nodes2)) 

022 (return nil)))))) 
023 

02* (defun equal' (desl des2) 

625 (cond < (eq desl des2) (list strue-node*) ) 

026 (t (insure-all-marked (list desl des2>) 

627 (if (eq (marker desl) (marker des2)) 

628 (list (eq-origin-node desl) (eq-or igin-node des2)))))) 
629 

636 (defun intern-equal' (desl des2) 

631 (if (eq desl des2) 

832 * true-node* 

833 (let ((expl (expansion desl)) 
83* (evp2 (expansion des2))> 

835 (if (and (not (atom expl)) 

836 (not (atom exp2)> 

837 (= (ncdrs expl) (ncdrs exp2))) 

638 (do ((desiqsl expl (cdr desigsl)) 

639 (desigs2 exp2 (cdr desigs2>> 
6*8 (nodes nil)) 

8*1 ((null desigsl) nodes) 

842 (let <(nodes2 (equal? (car desigsl) (car desigs2>>>> 

6*3 (if nodes2 

8** (setq nodes (nconc nodes nodes2)) 

6*5 (return nil)))))))) 

6*6 
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882 (defun * = (e*p) p>p:i 

883 (equality (designator expl) (designator exp2H) 
884 

885 (defun equality (r:j»sl na^e2> 

686 <)( (eq n^mol nam«2) 

887 Mru*»-node; 

888 (1st <*«wi-nodP (equality-node namel name?))) 

889 -~lcond ((not (evaluated? eq-node)) 
618 * _ (evaluated 1 eq-node) 

811 (addf (cons nano2 eq-node) (eq-links namel)) 

812 (add) (cons namel eq-node) (eq-links name2)) 

813 ( true-not icer eq-node '( I ink -check ' ,namel ',name2 ', eq-node)) 
81* (let <(»q-nodes (equal? namel name2))) 

815 ( t f eq-nodes 

616 (implies eq-nodes eq-node))))) 

617 eq-node) )> 

818 

819 (defun equality-node (desl des2) 

828 (dependency-node (if (eq 'less (scomp ( trans lat ion desl) (translat ion des2) ) ) 

821 ' ( = = .(translation desl) , (translat ion des2>) 

822 '( = = , (translation des2) , (translat ion desl) ))) ) 
823 

82* (defun I ink-check (desl des2 eq-node) 

825 (i f (marked^ desl) 

626 (i f (narked'' des2) 

827 (if (eq (marker desl) (marker des2>> 

628 (implies (list (eq-origin-node desl) 

829 (eq-origin-node des2>) 
838 eq-node) 

831 (i f (true -1 eq-node) 

832 (if (funcall Jhetter-pred* 

833 (canonical (marker desl)) 
83* (canonical (marker des2))) 

835 (mark des2 (marker desl) (eq-origin-node desl) eq-node) 

836 (mark desl (marker des2) (eq-origin-node des2) eq-node)))) 

837 . (if (true? eq-node) 

838 (mark des2 (marker desl) (eq-origin-node desl) eq-node))) 

839 (if (and (marked' des2) (true' eq-node)) 

8*8 (mark desl (marker des2) (eq-origin-node des2) eq-node)))) 

841 
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002 (defun ini t iflle-wrk«r (des) 

003 . (let ((ml (make-marker))) 
00* (setf (or iqin ml) des> 

005 (setf (canonical ml) des) 

006 (lei ((nodel (dependency-node '(canonical ,»1 , (translat ion d«s)))l 

007 (node2 (dependency-node '(active ,ml)))) 

008 f se I ( (canonical -node ml) nodel) 

009 "-.j- (set-truth nodel "(rue 'initiator) 

010 "(setf (active-node ml) node2> 

011 (set-truth node2 'true 'initiator) 

012 (mark-2 des ml strue-node*> 

013 ml))) 
01* 

015 (defun mark (de- ml eq-or igin-nodel eq-node) 

016 (let ( (eq-or iqm-node2 (equality-node des (origin ml)))) 

017 (implies (list eq-or igin-nodel eq-node) 

018 eq-or iqin-node2) 

019 (let ((can (canonical ml))) 

020 (it (funcall abetter -pred* des can) 

021 (make-canonical ml des)) 

022 (add! ' (.tcoincidencet ',can ',des) *think -queue*)) 

023 (mark-2 rt>5 ml eq-or igin-node2) ) ) 
024 

025 (defun make-canonical (ml des) 

026 (let ( (mot icers; ni I) 

027 (: removed- I is t? nil) 

028 (node (dopendoncy-node '(canonical ,ml , (translat ion des))))) 

029 (remove-2 (canon i ca I -node ml)) 

030 (setf (canonical ml) des) 

031 (sett (canonical-node ml) node) 

032 (set-2 node 'true 'premise 'marker) 

033 (removed-check) 

034 (run-not icers) ) ) 
035 

036 (defun mark-2 (des ml node) 

037 (let Um2 (marker des))) 

038 (l( (and m2 (not (eq ml m2>> (true' (active-node »2>)> 

039 (remove-truth (active-node m2)>)> 
0*0 (sett (marker des) ml) 

0*1 (sett (eq-or iqin-node des) node) 

0*2 (stimulate des) 

0*3 (mapc '(lambda (link) (link-check des (car link) (cdr link))) 

0** (eq-l inks des))) 

045 
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882 (defun stimulate (des) 

883 <M (or (null :name-queue*) 

884 (not (Junes I I *bet ter-pred* (car sname-queue*) des))) 
085 (add I des ^name-queue*) 

886 (insert dns *name-queuet) ) ) 

887 

888 (deftail insert (rte-s queue) 

089 (if (or Xnull (cdr queue)) 

818 (not (luncall tbet ter-pred* (cadr queue) des))) 

811 (addf des (cdr queue)) 

812 (insert des (cdr queue)))) 
813 

014 (deftail think () 

815 (cond (*■ think -queue? 

816 (let ((closure (car t th ink-queue*) ) > 

817 (relq > th ink -queue t (cdr *think -queue*) ) 

818 (|J think-trace (/Sprint ( I ist-translat ion closure))) 

819 (eval closure) 
828 (think))) 

821 (tname -queued 

822 (let ((des (car *name-queue*>) ) 

823 (setq ^name-queued (cdr tname-queue*)) 
024 (if (and (not (primitive 7 des)) 

825 (not (true' (sub -memo -node des)))) 

826 (substitute des)) 

827 (if (eq des (class-name des)) 

828 <think2 des)) 

829 (think))))) 
838 

831 (deftail think2 (des) 

832 (let ((closure (car (si imulate-l ist des)))) 

833 (cond (closure 

834 (delf closure (st imulate-l ist des)) 

835 (if think-trace (/fprint (list-translation closure))) 

836 (eval closure) 

837 (thmk2 des))))) 
838 

839 (defun print-val <evp> 

848 (let ((des (designator e*p)>) 

841 ( insure-marked des) 

842 (think) 

843 (/#prmt '((the value of) 

844 .(translation des) 

845 is 

846 .(translation (class-name des))) >)) 
847 
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082 (defun fund-co inc idence (desl des2) ~ ' "~ 

803 nil) 
084 

805 (defun fund-better (desl des2) 

806 (eq 'less (fund-comp desl des2>)) 
007 

008 (defun fund-comp (rtesl de--,2> 

089 (comp (suprr ioi-c lass 'unique' desl des2) 

010 'Ies5 

OH (comp (superior-class 'opaque' desl des2) 

012 'greater 

013 (num-comp (virt-weight desl) (virt-Meiqht des2>) 
0U Mess) 

015 'qreaterH 

016 

017 

018 

019 (defun unique' (des) 

020 (let ( (unp (uniqueness des))> 

021 (or (eq imp ' true) 

022 (and (not (eq unp 'false)) 

023 (cond ((numberp (eypansion des>) 
02* Isetf (uniqueness des) 'true) 

025 ,) 

026 ((pr imi t ive' des) 

027 (setf (uniqueness des) 'false) 

028 nil) 

029 ((a 1 1 -unique' (expansion des)) 
038 (setf (uniqueness des) 'true) 
831 t) 

®32 (t (setf (uniqueness des) 'false) 

033 nil)))))) 

834 

835 (defun all-unique' (desiqs) 

036 (or (null desiqs) 

037 (and (unique' (car desiqs)) 

038 (all-unique' (cdr desiqs))))) 
039 

840 (defun unique 1 (evp) 

041 (let ((des (designator exp>)> 

842 (if (uniqueness des) 

'*3 (break (unique! applied to an uniqueness determined object)) 

044 (setf (uniqueness des) 'true)))) 

045 
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888 
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817 

818 
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820 
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023 

024 

025 

026 

027 

028 

029 
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032 
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838 
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(def mac opaque' (des) 

(let ( (opq (opaqueness des))) 
(or (eq opq ' true) 

(and (not (eq opq ' false)) 
(cond ((primitive' des) 

(setf (opaqueness des) 'false) 
- nil) 

( (some-opaque' (expansion des)) 
(setf (opaqueness des) 'true) 
t) 
(t (setf (opaqueness des) 'false) 
nil)))))) 

(defun some-opaque" 1 (Hesiqs) 
(and des i qs 

(or (opaque' (car desiqs)) 

(some-opaque' (cdr desiqs))))) 

(defun opaque 1 (p.p) 

(let ((des (designator e*p))) 
t i I (opaqueness des) 

(break (opaque 1 applied to an opaqueness determined object)) 
(setf (opaqueness des) 'true)))) 



(defmac virt-weiqht (des) 
(let ((wgt (weight des))) 
(if wgt 
wgt 
(cond ((primitive' des) 

(setf (weiqht des) 1) 
1) 
(t (let ((uqt (sum-weight (expansion des)))) 
(setf (weight des) wgt) 
wgt)))))) 

(deftail sum-weight (def) 
(if (nul I del) 
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U (virt-weight (car def)) 
(sum-wetght (cdr def))))) 
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== EXPR 

»EXPAND EXPR 

»EXPHN0-2 EXPR 

ALL-UNIQUE? EXPR 

AUTO-SUB EXPR 

CANONICAL? EXPR 

CLASS-NAME EXPR 

DEFOPER DEFMACRO 

OESIGNRTOR EXPR .. 

DESIGNRTOR OEFTYPE 

DESIGNRT0R2 EXPR .. 

EQUAL' EXPR .. 

EQUALITY EXPR .. 

EQUALITY-NODE ... EXPR .. 

EQUIV-CLRSS EXPR .. 

EQUIV-CLRSS2 EXPR .. 

FUNO-BETTER EXPR .. 

FUNO-COINCIOENCE EXPR .. 



005 002 

oe: 03i 

002 039 
018 035 

003 020 

004 002 
003 014 
012 005 
002 002 
001 007 



002 
004 
005 
005 



005 
024 
005 
019 



007 004 
007 008 
018 005 
010 002 



FUND-COMP EXPR 010 808 

irWGE EXPR .... 083 032 

INITIATE-MARKER EXPR .... 006 002 

INSERT DEFTRIL 009 008 

INSURE-RIL-MRRKED EXPR 003 024 

INSURE -MRRKEO ... DEFHRCR0 883 811 

INTERN-CPNONICRL' EXPR .... 804 013 

INTERN-EQURL? ... EXPR 004 038 

LINK-CHECK EXPR 085 824 

MAKE-CANONICAL .. EXPR 886 825 

MARK EXPR 086 815 

HARK -2 EXPR .... 086 836 

MARKED' EXPR .... 003 882 

MARKER DEFSTRUCT 801 828 

NRME-INSPECT .... EXPR .... 882 023 

NRMES-IN1T EXPR .... 881 034 

NUMBER' DEFMACRO 001 025 



OPAQUE! EXPR . 

OPAQUE' DEFMOC 

OPERATOR-OEFS ... OEFflACRO 

PRIHIT1VE' DEFMACRO 

PRINT-VRL EXPR .. 

SOME -OPAQUE' EXPR .. 

STIMULATE EXPR .. 

SUBSTITUTE EXPR .. 

SUM-UE1GHT OEFTAIL 

THINK DEFTAIL 

THINK 2 DEFTRIL 

UNIQUE' EXPR .. 

UNIQUE' EXPR .. 

VIEU-ALL EXPR .. 

VIEH-CLRSS EXPR .. 

VIRT-UEIGHT OEFMRC 

UHRT-IS EXPR .. 



.. 811 


028 


.. Oil 


002 


812 


002 


001 


022 


.. 009 


039 


.. 811 


015 


.. 009 


002 


.. 808 


002 


Oil 


041 


009 


01* 


009 


031 


.. 010 


040 


.. 010 


019 


.. 007 


016 


.. 087 


019 


.. 011 


030 


.. 003 


018 
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Appendix Three 

The Algebraic Simplifier 

Algebraic simplification is done in the equality system in the same way 
as numerical computation. Evaluation functions are placed on algebraic operators 
which generate equalities between algebraic expressions and simplified forms of 
those expressions. Simplified versions of algebraic expressions are generated with 
a simplification function which operates on algebraic expressions and is 
completely separate from the equality system and the TMS. This section 
discusses the algorithms employed by that simplifier. 

There are some well known algorithms for handling the simplification of 
ratios of polynomials [Knuth 69]. These algorithms were implemented in a 
symbolic mathematical package developed at MIT [MACSYMA 77], and from 
there incorporated into a recent electrical circuit analysis system (SYN) developed 
by Johan de Kleer and Gerald Sussman [de Kleer & Sussman 78]. In the SYN 
system all expressions are converted to a canonical form which is a ratio of two 
polynomials. The variables are given global priorities and the polynomials in a 
canonical representation are polynomials in the highest priority variable occurring 
in them. The coefficients of such polynomials are canonical polynomials in the 
other variables. Synthetic division and greatest common divisor algorithms are 
used to remove common factors from the numerator and denominator of a ratio 
of two polynomials. This canonical form is truly canonical in that any two forms 
of a rational function will simplify to the same canonical expression. In this way 
simplification can be guaranteed to the extent that any expression can be placed 
in canonical form. 

A major problem with the symbolic algebra routines used in SYN is that 
computation of the greatest common divisor (gcd) of two polynomials in many 
variables is a very expensive operation. It was found in fact that this was the 
major limitation on the complexity of the circuits which could be analyzed in 
SYN. While recent developments in polynomial gcd techniques may improve this 
situation [Zippel 79], the simplification routines used in the equality system 
represent an attempt to avoid heavy reliance on gcd algorithms. 

Instead of reducing all expressions to a canonical polynomial form the 
algebraic simplifier used in the equality system is capable of maintaining factored 
forms of expressions. This allows factors to be removed from the numerator and 
denominator of a ratio by a simple search for common factors. It is not practical 
to factor all polynomials since this is much harder than computing geds, but once 
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an expression appears in factored form that factorization is not lost. This has the 
disadvantage that expressions are not placed in a truly canonical form, i.e. two 
expressions which arc really the same might be simplified differently if one was 
more factored than the other. 

In the equality system simplifier expressions are placed in pseudo- 
canonical form which is actually very similar to the canonical form used in SYN. 
Each expression in this pseudo-canonical form is a ratio of two products of 
polynomials. The variables have global priorities and polynomials are always 
polynomials in the highest priority variables appearing in them. The coefficients 
of these polynomials are again products of polynomials in lower priority variables 
and factors appearing in all the coefficients are factored out of the polynomials. 
Factors which are common to both the numerator and denominator of an 
expression arc removed from a ratio, but no attempt to find hidden common 
factors via a gcd algorithm has been incorporated. 

It would be very simple to add the ability to factor simple quadratics. It 
would also be possible to incorporate a gcd algorithm to guarantee complete 
removal of common factors from the coefficients and cancellations of all factors 
common to the numerator and denominator of a ratio. The polynomials to 
which the gcd would be applied would be considerably simpler in the system 
developed here than in a system which does not maintain factored forms. 

The simplification system uses an architecture developed by Gerald 
Sussman in which the basic simplification mechanism is a set of "symbolic 
operators" which compute the pseudo-canonical form for the sum, difference, 
product or ratio of expressions already in canonical form. A procedure for 
simplifying arbitrary expressions can be built which uses the symbolic operators to 
recursively simplifying subexpressions and then applies the symbolic operator 
corresponding to the top level operation. Products and ratios of expressions in 
pseudo-canonical form are very easy to handle since all that must be done is the 
merging of factors and the cancellation of factors common to the resulting 
numerator and denominator. Subtraction is trivially reduced to multiplication 
and addition, which is the most involved symbolic operation. 

Addition is done in two phases. First the expression is converted, in a 
straightforward way, to a ratio of a sum of two products and a product. The 
second step is the addition of the two products in the numerator. The first step 
in the addition of products is to remove common factors. Next, the highest 
priority variable of the two products is found, and each product is expanded into 
a polynomial in this variable. The coefficients of these polynomials are then 
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recursively added term by term. Finally factors which are common to all the 
coefficients of the resulting polynomial are removed a id the result of the original 
addition is a product of polynomials. 

If any product or polynomial can be simplified to zero this algorithm is 
guaranteed to do that simplification. This can be shown by observing that if any 
of the factors of a product simplify to zero then the product will also, and that 
all the factors are polynomials in some variable. A polynomial in a given variable 
will simplify to zero only if ail of its coefficients do, and these coefficients are 
either numeric or are simpler products of polynomials which, by an induction 
hypothesis, simplify to zero if such a simplification is possible. 

In the equality system each subexpression of an expression is a designator 
in its own right and can have properties attached to it. This can result in 
substantial space savings since all references to a particular expression point to 
the same data structure. This also makes it easy to memoize the simplification 
process which can result in significant time savings when an expression appears 
inside several different larger expressions or several times within the same 
expression. 

Another important point about the algebraic simplification system is that 
subexpressions which do not have top level algebraic operators are treated as 
variables. This allows the algebraic system to manipulate arbitrary functional 
expressions, which it has no knowledge of. For example there is no problem in 
dealing with (sin x) or (In (** x y)). These sub-expressions can be dealt with 
independently in the equality system, which is capable of substituting 
simplifications of them into algebraic expressions. 

The code for the algebraic simplification routines follow. 
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882 (defun operator (evp) 

883 (if (not (algebraic' exp)) 
804 nil 

085 (car evp))) 

886 

887 (defun operands (evp) 

888 (if (not laUqebraic' exp) > 
809 nil — 

818 (cdr »v P ))> 

011 

012 (defun make-prod (scale sumlist) 

013 (cons '* (cons scale sumlist))) 
014 

015 (defun numer-prod (evp) 

816 (cond ( (numbTp evp) (make-prod evp nil)) 

817 ((not (algebraic? evp)) (make-prod 1 (list exp))) 

818 (t (let ((op (operator exp))) 

819 (cond ( (eq op '// ) (numer-prod (cadr exp))) 

820 ( (eq op '*) exp) 

821 (t (make-prod 1 (list exp)))))))) 
822 

823 (defun denom-prod (evp) 

824 (cond ((not (algebraic 7 evp)) (make-prod 1 nil)) 

825 (t (let ((op (operator exp))) 

826 (cond ((eq op '// ) (numer-prod (caddr exp))) 

827 (t (make-prod 1 nil))))))) 
828 

829 (defun pnf_zero' (evp) 

830 ( = (car (operands (numer-prod exp))))) 
031 

032 (defun algebraic"' (e*p) 

033 (and (not (atom exp)) 

834 (not (numberp exp)) 

835 (memq (car exp) •(♦-*// )))) 
836 

837 (defun numerical' (e<p) 

838 (or (numberp evp) 

839 (let ((op (operator exp))) 

848 (cond ( (eq 't op) (null (cdr (operands exp)))) 

841 {(eq 'rat (operator evp>) t) 

842 ((eq '// (operator exp)) 

843 (and (numerical' (numer-prod exp)) 

844 (numerical' (denom-prod exp)))))) )) 
845 

846 (defun variable' (evp) 

847 . (not (or (algebraic' exp) 

848 (numerical' exp)))) 
849 

858 (defun sign (prod) 

851 (let <(n (cadr prod))) 

852 (make-prod (// n (abs n)) nil))) 
053 

854 (declare (special tvar-comp*)) 

855 

856 (setq *var-compt 'scomp) 

857 

858 (defun make-best-var (var old-var-comp) 

859 '(lambda (varl var2) 

868 (comp (superior-class '(lambda (var3) (eq var3 ',var)> 

861 varl var2) 

862 -less 

863 (, old-var-comp varl var2) 

864 'greater))) 
865 

066 (defun poly-comp (polyl poly2> 

867 (comp (funcall *var-comp* (poly-var polyl) (poly-var poly2)) 

868 'greater 

869 (scomp polyl poly2) 

870 Mess)) 
871 
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002 (defun *_«•_* (prodl prod?) "."'~ 

003 (let ((factorsl (operands prodl)) 
00* (factors? (operands prod?))) 

885 (iMKe-prod (t (csr fec-torsl) (car lactors2)) 

086 (comb-factors (cdr factorsl) (cdr factors2))))) 

007 

008 (defun comb-factors- (factorsl factors?) 

009 (cond (<r*ull factors]) factors?) 

010 ((null factors?) factorsl) 

811 (t (comp (poly-comp (car factorsl) (car factors2>) 

"12 (cons (car factors?) 

813 (comb-factors factorsl (cdr faetors2>>) 

01* (cons (car factorsl) 

815 (cons (car factors2) 

816 (comb-factors (cdr factorsl) (cdr factors?)))) 
017 (cons (car factorsl) 

818 (comb-factors (cdr factorsl) (actors2) )) )>) 
019 

020 (defun gcd_<_: (prortl prod2) 

021 (if (=0 (cadr prodl)) 

022 (make-prod 1 ni I ) 

023 (make-prod (qcd (cadr prodl) (cadr prod2>> 

82* (prod-intersection (cddr prodl) (cddr prod2))))> 

025 

02S (defun prod- inter-nct ion (factorsl factors2) 

027 (cond ((null factorsl) nil) 

028 ((null factors?) nil) 

829 (t (comp (poly-comp (car factorsl) (car faetors2>) 

830 (prod-intersection factorsl (cdr factors2>> 

831 (cons (car factorsl) 

832 (prod-intersection (cdr factorsl) (cdr f*ctors2)>) 

833 (prod-intersection (cdr factorsl) factors2) ) ) ) ) 
034 
835 
036 
037 (defun //_«_* (prodl prod2) 

838 (make-prod (// (cadr prodl) (cadr prod2>) 

839 (factor-deletion (cddr prodl) (cddr prod2>>>> 
040 
041 (defun factor-deletion (factorsl factors2) 

842 (cond ((null factorsl) nil) 

843 ((null factors2) factorsl) 

8** (t (comp (poly-comp (car factorsl) (car factors2>) 

8*5 (factor -deletion factorsl (cdr factors?)) 

8*6 (factor-deletion (cdr factorsl) (cdr tactors2)> 

847 (cons (car factorsl) 

8*8 (factor-deletion (cdr factorsl) «actors2))))>) 



; this division algorithm assumes that prodl is an explicit multiple of prod2. 
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(defun terms (poly) 

( i f (variable"* poly) 

(list (make-prod 1 (list poly))) 

(operands poly))) 

(defun poty-var (poly) 

(cond ( (mil I- po ly) nil) 

((variable' poly) poly) 
( lnu»?r ical' poiy) nil) 

(t (let ((first (cadr (operands (car (operands poly)))))) 
(cond ((variable? lirst) first) 

<t (break non-polynomial -sum)) ) )) ) ) 

(defun poly-mult (termsl terms? var) 

(cond ((null var) (break 'mil l-var-in-poly-mul t>) 
( (nul I termsl) nil) 
(t (let ((first (car termsl))) 

(poly-add (mapcar '(lambda (prod) (*_*_* first prod)) 
terms2) 
(poly-mult (cdr termsl) !erms2 var) 
var))))) 

(defun poly-add (termsl terms2 var) 

(let ((result (poly-add2 termsl terms2 var))) 
(cond ((pnf_zero ? (car result)) 
(cdr resul t>) 
(t result)))) 

(defun poly-add2 (termsl terms2 var) 

. (cond ((null var) (break 'nul l-var-in-poly-add) ) 
((null termsl) terms?) 
((null terms2) termsl) 

(t (let ((orderl (orrier2 (cdr (operands (car termsl))) var)) 
(order2 (order2 (cdr (operands (car terms2)>) var))) 
(cond ((= orderl order2) 

(cons (prodsum (numer-prod (car termsl)) 
(numer-prod (car terms2))) 
(poly-add (cdr termsl) (cdr terms2) var))) 
((> orderl order2) 

(cons (car termsl) (poly-add (cdr termsl) terms2 var))) 
(t (cons (car terms2) (poly-add termsl (cdr ter»s2> var)))))))) 

(defun order (poly) 

<erder2 (cdr (operands (car (operands poly)))) (poly-var poly))) 

(defun order2 (factors var) 

(if (eq (car factors) var) (1+ (order2 (cdr factors) var)) 8)) 



■ '*•*. *WW 3 W1W9 •««• • 
M2 (feftm efficients (poly) "~ " - 

Hi (cond ((variable' pely) (list (Mke-ere* * nit) 

••* (Mt«-»t-etf l nil))) 

Mf ((nurture eely) (list dm er ad »*!« Ml))) 

/""""N •••. <» (eoef2 (riwn (»»«ra«4« e«ly)> 

••' (»«lu-v*r eolu) 

••• t»»> 

■'-■■..•!• 

•4t (defun coef2 (terms var ini (-order) 

ill (if (mill (eras) 

•42 nil 

*** <*tt <(el (»rder2 (cdr ( eeeraw d a (c«r term) I) *•»»)) 

•1* (if (• et inM-erder) 

•4S (1*1 (l|-t«rH ( aa tr an ds (car teraalHI 

•J* («•»»• (Mete-erod tear **4*r*s> 

'■'•*' ***♦§ <«* .ft**er«*i IJ 

•*• «*»•« (eeV ter»** v«r #!♦ ta«*-erdtr>>)» 

•*• <e«nt • (eeef2 teres vtr «♦ inltHtrdl 

•2t 

•24 (defun coef3 lltrtt var> 

•22 (if («q (car tere*) var) 

•29 (coe(3 (cdr Mw) ver) 

•24 teres)) 



^ 
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W2 (datun -_pnf_l (e«p) 
M3 (lat (<n (miwar-prod axp>>) 

••4 (Mat '// 

_^^ M5 (Mke-prad <* -1 (cadr n» (cddr n>) 

r \ ••• (dana*-pre<l axp>>>> 

•17 

Mt (dalun -_pnf (argsr 
•ft (it (cdr-arqs) 

•1* («jmf~<Mtt (ear argt) 

•it (-j»f«t_l (eaa> «-«•>>>> 

•i2 <-_P»»«_l (car ar«>>>> 

•13 

•14 (dafim \/t («xp) 

•15 <l>«t *// (o>no»-pra4 axp) (mmtr-arod exp>>> 

•» 

•17 (da fun //_pnl («rqs) 

•It (»_pn( (list (car aras> (1// (cadr aras>»>> 

•19 

•*• (da tun *_pnt (arg«> 

•21 (cond (I* (Itngth ar«*> 1) (ear arja)) 

•22 (t (tiaprad (ear aros) (*_pnt (eaV arvsH)))) 

•23 

•2* (datun siaprod (evpl axp2> 

•25 (lat ((top-prod (*_*_» (maiar-pred axpll 

•28 (nuaar-pred exp2)M 

•27 (botto*-prad (*_*_* (dana»-prad axpl) 

•2» ( dan e * prad axp2>»> 

•2f (tat ((coanon-prad (fed.*-* tap-prad a*ttaa). ar a d» > 

•3« (list '// 

•31 <//_♦_* tap-pred cowa n pi a O 

•32 (//.*_* battaacpred co— aw prad»»>> 

•34 (datun ♦_**»( (args) 

•35 (cond ((> (lanajh area) D (car arcs)) 

•3» (t (siapatm (ear ar«<) (♦ J «» (cdr araaUHH 

•37 
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882 (defun simpsum (expl exp2> 

883 (let <(»opl (t_t_* (numer-prod expl) (denom-prod exp2))) 

884 Hop2 <«_f_* (numer-prod exp2) (denom-prod expl)))) 
085 (let ((top-prod (prod^uni topi top2)) 

886 (bottom-prod (t_s_* (denom-prod expl) (denom-prod exp2>))) 

887 (let ((common-prod (qcd_f_* top-prod bottom-prod) )) 

888 . (if (not (pnfzero? top-prod)) 

889 __ (ii S t V/ 

818 - '//_*_* top-prod common-prod) 

811 '//_*_* bottom-prod common-prod)) 

812 (if (pnf_zero' bottom-prod) 

813 (break indeterminate-value) 
81* (list '// top-prod 1))))))) 
815 

816 (defun prodsum (prodl prod2) 

817 (let ((corn-fact (corn-factor (list prodl prod2)>>> 

818 (let ((rest-prod (con-sum (//_«•_* prodl corn-fact) 

819 <//_*!* prod2 corn-fact)))) 
828 (if (s (car (operands rest-prod)) 8) 

821 (make-prod 8 ni I) 

822 (i t_* corn-fact rest-prod))))) 
823 

824 (delun con-sum (terml term?) 

825 • (if (and (numerical 7 terml) (numerical' term2>) 

826 (make-prod (♦ (car (operands terml)) (car (operands term2))) nil) 

827 (let ((varl (poly-var (cadr (operands terml)))) 

828 (varZ (poly-var (cadr (operands term2))))) 

829 (cond ((eq varl var2) 

838 (factor (poly-add (expand terml) (expand tern>2) varl))) 

831 ((or (null var2) (and varl (eq 'less 

832 (funcall *var-comp* varl var2>>>) 

833 (factor (poly-add (expand terml) (list term2) varl))) 

834 (t (factor (poly-add (expand term2) (list terml) var2)>>))>) 
835 

836 (defun expand (prod) 

837 (let ((opers (operands prod))) 

838 (if (null (cdr opers)) 

839 (I 1st prod) 

848 (evpand2 (car opers) (cdr opers))))) 
841 

842 (defun expanH2 (scale factors) 

843 (cond ((null factors) (list (make-prod scale ni I)) ) 

844 ((eq (poly-var (car factors)) (poly-var (cadr factors))) 

845 t (poly-mult (terms (car factors)) 

8*6 (expand2 scale (cdr factors)) 

847 (poly-var (car factors)))) 

'848 (t (let ((rest (make-prod scale (cdr factors)))) 

849 (mapcar '(lambda (prod) (*_*_* prod rest)) 
858 (terms (car factors))))))) 

851 

852 (defun factor (terms) 

853 (cond ((null terms) (make-prod 8 nil)) 

854 ((null (cdr terms)) (car terms)) 

855 (t (let ((common (■*_*_* (sign (car terms)) (corn-factor terms)))) 

856 '*_*_* common 

857 (make-prod 1 (list (cons '♦ (mapcar '(lambda (prod) 

858 <//_*_» prod common)) 

859 terms)))))))))" 
868 

861 (defun corn-factor (sumlist) 

862 (com-factor2 (car sumlist) (cdr sumlist))) 
863 

864 (defun com-factor2 (factor sumlist) 

865 (cond ((null sumlist) factor) 

866 (t (com-factor2 (gcd_*_» factor (car sumlist)) (cdr sumlist))))) 
067 
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882 jsolve returns a dotted pair. 

083 j the car ol the pair is a list of sums uhoes product must not be zero. 

88* ;the cdr ot the pair is a list ol roots ol the equaltion, they are not In canonical 

885 ; form however and must be sir.pl ified before being passed as arguments to symop functions. 

886 

887 (delun solve (evp var) 

808 _ (let* ( (;v,ir--comp ; <make-b?st-var var *var-comp*)> 

889 (e»p2 (simplification e*p>>> 

810 (if lor (pnlzero 7 (numer-prod e*p2)) 

011 (pn)_z«ro' (denom-prod exp2))) 

012 nil 

813 (let (<(non-z«>ro . roots) 

"1* (solved (cdr (operands (numer-prod exp2))> var))) 

815 (cons (append (mapcar 'alg-trans (cdr (operands (denom-prod exp2)>)) 

016 non-zero) 

017 roots))))) 
018 

019 (delun solve? (factors var) 

020 (if factors 

021 (let (((non-zero . roots) (solve2 (cdr factors) var)) 

022 (sum (car factors))) 

823 (if (eq var (poly-var sum)) 

02* (let* ((coels (coelicients sum))) 

&25 (cons (append (mapcar 'alg-trans (cdr (operands (car (last coefs))))) 

826 non-zero) 

827 (cons (if (cddr coefs) 

828 '(a-root-of ,§(mapcar 'alg-trans coefs)) 

829 (alg-trans (//_pnf (list (-_pnf_l (car coefs)) 

830 (cadr coefs))))) 
031 roots))) 

832 (cons (cons (alg-trans sum) non-zero) 

833 roots))))) 
834 
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002 (defun variables l»vp) 

083 (merge (var iables-2 exp) nil)) 

004 

085 (defun variables-? (exp) 

006 (cond ((alqebraic? exp) 

007 (mapcan 'variables-2 (operands exp))) 

008 ( (mimsr ical ' exp) nil) 

009 <t_llist exp)))) 
010 

011 (defun simplification (evp) 

012 (if (not (algebraic' exp)) 

013 evp 

81* (let ((args (mapcar 'simplification (operands exp))) 

015 (op (operator exp))) 

016 (cond ( (eq op '♦) (+_pnf args)) 

817 (<eq op '-) (-_pnf args)) 

818 ((eq op '*> (*_pnf args)) 

819 ((eq op '// ) "(//_pnf args)))))) 
020 

021 ;the next function makes expressions more readible by removinq factors of 1 etc. 

022 

023 (defun alg-trans (evp) 

824 (cond ((not (algebraic' exp)) 

025 evp) 

826 (t (let ((op (operator exp)) 

827 (opers (mapcar 'alg-trans (operands exp)))) 

828 (cond ((eq op '// ) 

829 (cond ((equal (cadr opers) 1) (car opers)) 

830 ((equal (cadr opers) -1) 

831 (cond ( (numberp (car opers)) 

832 (* -l (car opers))) 

833 ((not (algebraic' (car opers))) 

834 '<* -1 , (car opers))) 

835 ( (eq (operator (car opers)) '*) 

836 •(* , ($ -i (cadar opers)) , (caddar opers))) 

837 (t '<* -1 ,(car opers))))) 

838 (1 (cons '// opers)))) 

839 ((eq op '4) (cons '♦ opers)) 

840 ((eq op '*) 

8*1 (cond ((equal (car opers) 0) 0) 

8*2 ((equal (car opers) 1) 

8*3 , (cond ((null (cdr opers)) 

044 1, 

8*5 ((null (cddr opers)) 

8*6 (cadr opers)) 

8*7 • (t (cons '* (cdr opers))))) 

8*8 ((nul I (cdr opers)) 

8*9 (car opers) > 

858 (t (cons '* opers))))))))) 

051 

852 (defun simpt (exp) 

853 (alg-trans (simplification exp))) 
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*VAR-C0f1P* SETQ 001 056 

*_*_* EXPR oo: 002 

*_pnf expr cos oro 

♦_PNF EXPR 006 034 

-_PNF EXPR 006 008 

-_PNF_1 EXPR _§06 002 

/_*_# EXPR 042 037 

/_PNF EXPR 006 017 

l7 EXPR 006 OH 

AIG-TRANS EXPR 009 023 

ALGEBRAIC EXPR 001 032 

C0EF2 EXPR 004 010 

C0EF3 EXPR 004 021 

COEFICIENTS EXPR 004 002 

COrt-FBCTOR EXPR 007 061 

C0t1-FflCT0R2 EXPR 007 064 

COHB-FflCTORS EXPR 002 008 



CON-SUM EXPR 

DENOn-PROO EXPR 

EXPRNO EXPR 

EXPRN02 EXPR 

FACTOR EXPR 

FRCTOR-DELETION EXPR 

GCO_:_: EXPR 

MRI'E-flEST-VRR ... EXPR 

nni e-ppoo expr 

NUMFJR-PROO EXPR 

NUMERICAL' EXPR 

OPERANDS EXPR 

OPERATOR EXPR 

ORDER EXPR 

0RDER2 EXPR 

PNF_ZERO' EXPR 

POLY-ADO EXPR 



007 


824 


001 


023 


007 


036 


007 


042 


007 


0S2 


002 


041 


002 


820 


001 


858 


001 


012 


001 


015 


001 


037 


001 


007 


001 


882 


003 


844 


003 


847 


001 


829 


803 


82* 



P0LY-RD02 EXPR 003 030 

POLY-COrtP EXPR 001 066 

POLY-MULT EXPR 003 015 

POLY-VRR EXPR 003 007 

PROD-INTERSECTION EXPR 002 026 

PROOSUri EXPR 007 016 

SIGN EXPR 001 050 

SIMPLIFICATION .. EXPR 009 Oil 

SIMPROO EXPR 806 024 

SIMPSUfl EXPR 007 002 

SIMPT EXPR 009 052 

SOLVE EXPR 008 007 

S0LVE2 EXPR 808 019 

TERMS EXPR 083 002 

VARIABLE? EXPR 001 046 

VARIABLES EXPR 009 002 

VARIABLES-2 EXPR 809 005 
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802 (declare (special tzero* plunic-we ight-count)) 

883 

884 (defun algsys-init () 

885 (names-init ' a lysys-bet ter 'algsys-coin) 

886 (mapc 'unique 1 '<♦-//*)> 

887 (setq ;;?roi (designator 0>) 
008 (name- type "_p I untc ) 

089 (setq plunk -weight -count D) 

818 

811 (putprop '♦ '*_pnf 'symop) 

612 (putprop '- *-_pnf 'symop) 

613 (putprop ' ■» 't_pn( 'symop) 
81* (putprop '// V/.pnf "symop) 

815 (addf 'create-s imp I 1 f icat ion (operator-dels '♦)) 

818 (addf 'create-simpi i f icat ion (operator -defs '-)) 

817 (addf 'create-s imp I j f teat ton (operator-defs '*)) 

818 (addf 'create-simpi i f icat ion (operator-defs ' // )) 
019 

028 (defun create-s impl i f icat ion (des) 

021 (let ((simp (des-simpl i f icat ion des) )) 

822 (let ((desT (designator2 (alg-trans simp)))) 

823 (associate 'simplification simp (memoizat ions des)) 
024 (set-truth (equality des des2) 'true 'alqebra)))) 
825 

026 (defun des-s impl i f icat ion (des) 

827 (let < (exp (expansion des))) 

828 (cond ((atom exp) exp) 

829 ((not (memq (expansion (car exp)) '(+-// *)>) 
838 (I ist-translation des)) 

031 (t (let ((simp (assq 'simplification (memo i rat ions des)))) 

832 ( i f s imp 

833 (cdr simp) 

034 (let ((result (funcall (get (expansion (car exp)) 'symop) 

035 (mapcar 'des-simpl I f icat ion (cdr exp))))) 

036 (addf (cons 'simpl i f icat ion resul t> 

037 (memo izat ions des)) 
838 result))))))) 

839 

840 (defun c== (x y> 

0*1 (let* ((z-exp (alg-trans (des-simpl i f icat ion (designator '(- ,x ,y))))) 

0*2 (c-node (equality (designator2 z-exp) szero*))) 

0*3 (mapc '(lamhda (var) (solve-constraint z-exp var c-node)) 

0** (variables z-exp)) 

045 c-node)) 

846 

847 (defun solve-constraint (exp var c-node) 

8*8 (let ((t-var (designator var)) 

0*9 ((non-zero . roots) (solve exp var))) 

058 (if (and (car roots) (null (cdr roots))) 

851 (if (null non-zero) 

052 (implies (list c-node) 

053 (equality t-var (designator2 (car roots)))) 
05* (let* ((non-z (if (cdr non-zero) 

055 (designator2 '(* ,€non-zero)) 

056 (designator2 (car non-zero)))) 

057 <eq-node (equality non-z tzero*))) 

058 (default eq-node 'false) 

059 (implies (list (t-not eq-node) c-node) 

868 (equality t-var (designator (car roots)) )>))>) ) 
861 



001 fiPPNOX 4 18/21/79 Page 2 

802 (defun a Igsys-bet ter (desl des2) 

883 (eq 'less (algsys-comp desl des2>)> 

884 

805 (defun algsys-comp (desl dss2) 

886 (comp (super ior-c lass 'unique' desl des2) 

087 'less 

888 (comp .(super ior-c lass 'opaque' desl des2) 

089 'greater 

018 - (comp (num-comp (virt-plunk-ue igM desl) 

011 (virt-plunk-weigbt des2)> 

012 'greater 

013 (num-comp (virt-weight desl) (virt-weight des2)) 

014 'less) 

015 Mess) 

016 'qreater)) 
017 
818 
819 
028 
821 bpw 

823 (set! (plunk-weight des) 0) 

824 0' 

825 (t 
026 

027 hpu)))))> 

028 

829 (defun biggest-plunk -weight (desigs) 
038 (if (nuli (cdr desigs)) 

031 (vir t-plunk-we ight (car desigs)) 

032 (let (<bpul (biggest -plunk -weight (cdr desigs))) 

033 fhnu? t \r ■ r \ _« I unlr _ne ■ nk * f^ar> ftBeinellll 

034 
835 
836 
037 (defun plunk' (evp) 

MB / :_• « /_ i *. : 



(defun virt-plunk -weight (des) 
(let ((bpw (plunk -weight des))) 
( i f bpw 
»i bpw 

022 (cond ((primitive' des) 

""'' , (sell (plunk -weighl 

0) 

t (let ((bpw (biggest-plunk -weight (expansion des)))) 
(setf (plunk-weight des) bpw) 
bpw)))))) 

plunk-weight (desigs) 

'r desigs)) 

nk-weight (car desigs)) 
et (<bpwl (biggest -plunk -weight (cdr desigs))) 

<bpw2 (vir t-plunk-we ight (car desigs)))) 
(if (not (and bpwl bpw2>) (break (null plunk weight))) 
(if (> bpwl bpw2) bpwl bp«2>>>) 



040 

041 (se 

042 (in 

043 (se 
044 
845 



038 (print * (plunk inq , evp>) 

039 (let ((desl (designator exp)) 
""* <des2 (designator (intern (gename 'plunk))))) 

itf (plunk-weight des2) plunk-weigfrt-count) 

icremont plunk -we ight-count) 

it-truth (equality desl des2) 'true 'definition))) 



eel RPPNDX * 18/21/79 Page 3 

082 (defun alqsys-coin (desl des?) 

083 (if (eq desl dfs2) (break (eq dess in algsys-coin))) 
884 (if (and (eq desl (image desl)) 

005 (eq des2 (image des2))) 

086 (let ((t-desl (translation desl)) 

007 (t-des2 (translation des2)>) 

088 (cond ((and (numberp t-desl) (numberp t-des2)) 

009 -- (set-truth (equality desl des2) 'false 'algebra)) 

010 " ((and (numerical' t-desl) 

011 (numerical' t-des2>> 

012 (let <(diff (des-simpl if ication (designator '(- ,desl ,des2))))) 

013 (il (not (pnf_zero' diH)) 

01* (set-truth (equality desl des2> 'false 'algebra)))) 

015 ((or (> (virt-plunk-weight desl) 8) 

016 (> (virt-plunk-weight des2) 8)) 

017 (handle-plunk (designator 

018 . (alg-trans 

019 (des-s impl i f icat ion 

020 (designator '(- ,desl ,des2))))) 

021 (equality desl des2)))>>>> 
022 

023 (defun handle-plunk (des c-node) 

024 (let ((plunk (base-plunk des (virt-plunk -weight des)))) 

025 (if plunk (sol ve-constraint (translation des) 

026 ( trans lat ion plunk) 

027 c-node)))) 
028 

029 (defun base-plunk (des ueight) 

030 (if (primitive' des) 

831 (if (= weight (virt-plunk -ueight des)) 

032 des) 

033 (base-plunk2 (expansion des) weight))) 
034 

035 (defun base-plunk2 (desigs ueight) 

036 (if desigs 

037 (let ( (bp (base-plunk (car desigs) weight))) 

038 (if bp bp (base-plunk2 (cdr desigs) weight))))) 
039 

040 
041 
042 

043 (defun solve-for (evp) 

044 (think) 

045 (let ((des (desiqnalor evp))) 

046 ' (let ((can (class-name des) )) 

047 (if (and (not (opaque' can)) < = (virt-plunk-weight can) 8>) 
648 (translation can) 

049 (let ((desigs (equiv-class des))) 

850 (do ((rest desigs (cdr rest))) 

051 ((null rest) '(it does not seem possible to do so)) 

052 (if (and (opaque' (car rest)) 

053 (algebraic' (translation (car rest)))) 
854 (progn (mapc 'plunk 1 (opaque-vars (car rest))) 

055 (think) 

056 (return (translation (class-name des))))))))))) 
057 

058 (defun opaque-vars (des) 

059 (fmapcar '(lambda (var) 

060 (let ((des (designator var))) 

061 (if (opaque' des) des))) 

062 (variables (translat ion des)))) 
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Appendix Five 

Utility Procedures 

Basic utility functions and macros are described here which are used by 
the code in appendix two. The LISP definition of most of these utilities is given 
here. However only a brief english description of some of the basic utilities. 
Most of the basic concepts behind the utilities described here were developed by 
various people other than the author and many of them are documented in the 
LISP MACHINE MANUAL [Wienreb & Moon 78]. 

If 

(if a id c) is equivalent to: (cond (a b) (t c)). (if a b) is equivalent 
to: (cond (a b)). 

Backquote 

The backquote feature provides a form of quote which replaces items 
preceded by a comma with their value. The following are some examples of the 
use of backquote: 

'(foo a , (+ 1 2)) evaluates to: (foo a 3) 

Mfoo .(list 'a *b) (list 'a 'b)) evaluates to: (foo (a b) (list 'a *b)) 

Items in the interior of backquoted expressions which are preceded by ,e 
have there values exploded into the top level list structure. An example of the 
use of this feature is as follows: 

*(foo .©(list 'a *b) .(list 'a *b) (list 'a *b)) 

evaluates to: 
(foo a b (a b) (I ist 'a 'b)> 

Defmacro 

This form is used to define macros. A macro definition has a similar 
syntax to a function definition. When a form whose car is a macro is evaluated 
the macro definition is used to generate a new form whose value is the value 
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returned for the original form. The arguments to the macro are bound to the 
forms in the argument positions rather than their values as is done for functions. 
An example of a macro definition is given below: 

(defmacro first-part (x) 
'(caar ,x)) 

Using this definition (first-part a) macro expands to: (caar a) and so 
(fist-part a) has the same value as: (caar a). A macro is often used instead 
of a trivial function definition because it is expanded within the compiler and 
results in more efficient compiled code. 

It is sometimes convenient to allow the bound variable list of a macro to 
be an arbitrary list structure rather than a simple list. In this case atoms in the 
bound variable list (or bound variable pattern, since it need not be a simple list) 
are bound to corresponding parts of the expression using the macro. For example 
the new MACLISP form of do could have been defined as a macro along the 
following lines: 

(defmacro do (variable-bindings (end-test . end-body) . do-body) 

, The bound variable list may also be a single atom, in which case that 

atom is bound to the entire list of "arguments" to the macro. 

Defmac 

defmac is identical to defun with the exception that a macro is created 
which the compiler can use to open code the function during the compilation of 
other functions. This is used purely for reasons of efficiency. The open coding 
is useful in getting the compiler (and other optimization macros such as def tai i) 
to perform optimizations which would not otherwise be done. No function 
defined via defmac can be recursive however since this would lead to infinite 
expansion during open coding. 

Let 

The let feature allows structured lambda binding. An example follows: 
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(let {(a 1) 
(b 2)) 
(+ a b)) 

i s equivalent to: 

((lambda (a bl It a b)) 1 2) 

The let macro allows the the bindee of a binding pair to be an 
arbitrary list structure whose parts are bound to the corresponding parts of the 
value being bound. This is convenient for dealing with functions which 
conceptually return more than one value. 

Setf 

The setf macro gives a general method for side effecting data 
structures. The following equivalences give some examples of its use: 

(setf a b) is equivalent to: (setq a b) 

(setf (get a b) c) (putprop a c b) 

(setf (car a) b) (rplaca a b) 

(setf (cdr a) b) (rplacd a b) 

(setf (cond (a b) (c d)) e) (cond (a (setf b e) ) (c (setf d e) ) ) 

Defsidmac 

defsidmac is just like defmacro except that it is used to define macros 
which side effect their last argument and treats that argument position specially. 
Specifically it defines a macro which will embed the side effect in conditionals as 
does setfaddf. 

(defsidmac addf (x list) 

'(setf .list (cons ,x . list))) 

(addf x b) is equivalent to: (setf b (cons a b) ) 
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but 

(addf x (if a b c)) is equivalent to: (if a (addf x b) (addf x c) ) 

While it may seem obscure to write code which side effects conditional 
expressions, the ability to do so can be important when macros expand to 
conditionals. In such situations it is sometimes convenient to be able to side 
effect applications of these macros. 

Defstruct 

The defstruct feature is used to define a type of structured object. A 
defstruct definition creates a set of macros. One of these macros is used to 
create objects of the defined type. The others are used to access the parts of 
that object. Consider the following example: 

(defstruct (ship) x-pos y-pos (mass 280)) 

This defines four macros: make-ship, x-pos, y-pos, and mass. The 
make-ship macro creates a ship with its mass set to a default value of 200. The 
following dialogue illustrates a use of these macros: 

(SETQ HERO (MAKE-SHIP)) 
(ni I ni I 2801 

(MASS HERO) 
200 

(SETF (X-POS HERO) 10) 
10 

(X-POS HERO) 
10 

Deftype 

deftype is just like defstruct except that it defines another macro 
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which determines if any given object is of the defined type. In the above 
example for clef struct if deftype had been used instead then a ship? macro 
would have been created which could determine if any given object was created 

via make-ship. 

Deftail 

This feature allows automatic tail recursion optimization. The way this 
is invoked is to use deftail instead of defun in a function definition. This 
feature actually docs more than simple tail recursion optimization in that simple 
accumulations (functions which generate sums, products, or lists recursively) are 
also converted to iterative forms. The reader need not pay any attention to this 
feature and may assume that all functions defined via deftail are actually 
defined via defun in the standard fashion. 

Defarb 

defarb is identical to defun except that it allows the bound variable list 
to be an arbitrary list expression. The atoms in this expression are bound to the 
corresponding parts of the list of values to which the defined is applied. The 
most common use of defarb is to have the bound variable pattern be a single 
atom in which case that atom is bound to the list of arguments to the function. 
A function so defined can take an arbitrary number of arguments. 

The code for other utility functions and macros follows 
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(setf (cond (a b) (c d))) 

expand to: 
(cond (a (set* b e)> (c (setf d «))) 



;defsidmac defines cond- imbedding side effect macros I ike setf 
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002 (defmacro setf (form value) 

003 (let ((form (macroexpand form))) 

004 (cond ((symbolp form) '(setq .form .value)) 

005 ((numberp form) (break (attempt to set a number))) 

006 <<ea (car form) 'cxr) 

007 '(rplacx , (cadr form) , (caddr form) .value)) 

008 _ ((eq (car form) 'car) 

00 9 - Mrplaca , (cadr form) .value)) 

1 (<eq (car form) 'cdr) 

01 * '(rplacd , (cadr form) .value)) 

012 ( (eq (car form) 'get) 

013 "(putprop , (cadr form) .value , (caddr form))) 

14 ((eq (car form) 'cond) 

15 (cond-setf form value)) 

016 <t (break (unknown operator In setf)))))) 
017 
018 
019 
020 
021 

022 (defun cond-setf (form result) 

023 (cond- imbed '(lambda (form) '(setf , form , '.resul t)> 

024 form)) 
025 

026 (defun cond- imbed (side-effect cond-form) 

027 '(cond ,e(mapcar '(lambda (clause) 
°28 (if (or (cddr clause) (null (cdr clause))) 

029 (break (illegal cond clause for cond side effect)) 

030 '(.(car clause) 

031 .(funcall side-effect (cadr clause))))) 

032 (cdr cond-form)))) 
033 
034 
035 

036 (defmacro defsidmac (effect vars . body) 

037 (let ((rvars (reverse vars))) 

038 (let (<ref (car rvars)) 

039 (params (nreverse (cdr rvars)))) 

040 '(defmacro .effect ,vars 

41 (let ((ref2 (macroexpand ,ref))) 

42 (if (and (not (atomref2)> 

043 (eq (car ref2) 'cond)) 

44 (cond-imbed '(lambda (ref) '(,',*, effect 

045 ,e\(list .eparams) 
04S ,reO) 

047 ref2 ) 

048 (progn ,6body) >)) )) ) 
049 

050 (defsidmac increment (x) 

051 '(setf ,x (1+ ,x))) 
052 

053 (defsidmac addf (val list-ref) 
°54 '(setf , list-ref (cons ,val , list-ref))) 

055 

056 (defsidmac delf (val list-ref) 

057 '(setf , list-ref (delete ,val , list-ref))) 
058 

059 (defsidmac associate (x y assoc-ref) 

060 (let ((argl (namcopy 'x>) 

061 (arg2 (namcopy 'y)) 
G2 (alist (namcopy 'alist)) 

063 (as (namcopy 'as))) 

064 '(let ((.argl ,x) 

065 (,arg2 ,y) 

66 (, alist , assoc-ref)) 

67 (let ((.as (assoc ,argl , alist))) 

068 (if ,as 

69 (rplacd ,as ,arg2) 

070 <setf , assoc-ref (cons (cons ,argl ,arg2) ,al ist)) >)))) 
871 
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(delmacro hash-array (array-name dim) 

'(array , (cadr array-name) t ,dim)) 

(defun hash (item range) 

(remainder (abs (ranhash item)) range)) 

(deMail ranhash (item) 
(cond ((null i tem)"-0) 

((symbolp item) (symbol-hash item)) 
( (numberp item) (fix item)) 
((designator? item) (desig-hash item)) 
((node 7 item) (node-hash item)) 
( ( I istp i tem) 
(logxor (ranhash (car item)) 

(ranhash (cdr item)))) 
(t 0))) 

(defun node-hash (node) 

(let ((val (node-hash-val node))) 
(cond (val val) 

(t (setq val (random)) 

(set! (node-hash-val node) val) 
va I ) ) ) ) 

(defun desig-hash (desig) 

(let ((val (desig-hash-val desig))) 
(cond (val val) 

(t (setq val (random)) 

(setf (desig-hash-val desig) val) 
va I ) ) ) ) 

(defun symbol -hash (symb) 

(let ((val (get symb 'ranhash))) 
(cond (val val) 

(t (setq val (random)) 

(putprop symb val 'ranhash) 
va I ) ) ) ) 

(defun table-asscr (item table) 

(and (eq (typep table) 'symbol) (setq table (get table 'array))) 
(let* ((index (hash item (cadr (arraydims table)))) 
(bucket (arraycall t table index)) 
(result (assoc i tern bucket) )) 
(cond (result result) 
(t 
(setq result (cons item nil)) 

(store (arraycall t table index) (cons result bucket)) 
result)))) 
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023 

024 
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027 

028 

029 

030 

031 

032 

033 

034 
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(declare (special *query-stack«)) 

(defun push-query (query) 

(w-push query iquery-stack*) 
(vieu-query (w-top *query-stack*) ) ) 

(defun pop-query ()_ 

(w-pop *query-stack*> 

(vieM-query (w-top squery-stacK*) ) ) 

(defun view-query (query) 
(cons (car query) 

(mapcar '(lambda (n ass) (cons n (car ass))> 

(integers-betMeen 1 (length (cdr query))) 
(cdr query)))) 

(defun answer (n) 

(cdr (nth (1- n) (cdr (w-top *query-stack*>))>> 

(defun make-wrap-stack (depth) 

(let ((last (cons nil (cons nil nil)))) 
(let ((first (mws2 last depth))) 
(setf (car (cdr first)) last) 
(setf (cdr (cdr last)) first) 
first))) 

(defun mws2 (last depth) 
(if (= depth 1) 
last 
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(let 



((second (tnws2 last (1- depth))) 

(first (cons nil (cons nil nil)))) 
(setf (cdr (cdr first)) second) 
(setf (car (cdr second)) first) 
first))) 



(defsidmac w-push (item wstack) 

'(let <(wstack2 (cadr .wstack))) 
(setf (car wstack2) , i tern) 
(setf , wstack wstack2))> 

(defsidmac w-pop (wstack) 

'(setf .wstack (cddr , wstack))) 

(def macro w-top (wstack) 
'(car .wstack)) 



(progn (setq *query-stack* (make-wrap-stack 38)) 1) jprogn prevents inf inite print i 



ng 
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002 j the following macro is useful! in dealing with partial orders in Mhlch three conditional 

003 ; branches exist depending on the relation of the two items. 
004 

005 

006 (de (macro comp (comparison greater-case unrel-case less-case) 

007 (let ( (comp-var (namcopy 'comp-var))) 

008 '(let ((, comp-var .comparison)) 

009 (cond T\eq , comp-var 'greater) , greater-case) 

010 ((eq , comp-var Mess) , less-case) 

011 (t ,unrel-case))))) 
012 

013 (defmac num-comp (nl n2> 

014 (cond ((> nl n2> 'greater) 

015 <U nl n2> ' less) 

016 (t 'equal))) 
017 

018 (defmac alpha-comp (al a2) 

019 (cond ((alphalessp al a2) Mess) 

020 ((equal al a2> 'unrelated) 

021 (t 'greater))) 
022 

023 (defun scomp (expl exp2) 

024 (comp (num-comp (ncdrs expl) (ncdrs exp2)) 

025 'greater 

026 (lexical -comp expl exp2) 

027 Mess)) 
028 

029 (deftai I ncdrs (exp) 

030 (if (atom exp) (U (ncdrs (cdr exp))))) 
031 

032 (defun lexical -comp (expl exp2) ; they must have the same ncdrs 

033 (cond ((or (numberp expl) (numberp exp2>) 

034 (cond ((not (numberp exp2)) Mess) 

035 ((not (numberp expl)) 'greater) 

036 (t (num-comp expl exp2>))) 

037 ((atom expl) ;if one is then they both art 

038 (alpha-comp expl exp2)> 

039 (t (comp (scomp (car expl) (car exp2)> 

048 'greater 

0*1 (lexical -comp (cdr expl) (cdr exp2)) 

0*2 Mess)))) 

043 

844 (defmacro superior-class (pred x y) 

045 Mif (funcallm ,pred ,x> 

046 (if (funcallm , pred ,y> 

047 'unrelated 
848 'greater) 

049 (if (funcallm ,pred ,y) 

850 • less 

851 'unrelated))) 



Symbol Table for: 
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ROOF DEFS IOMRC 002 853 

RLPHR-COflP DEFMRC . . 085 018 

^NSUER EXPR 004 018 

iSSOCIRTE DEFS IOMflC 082 859 

COMP DEFMRCRO 005"0O6 

CONO-iriBEO EXPR "TO2 026 

CONO-SETF EXPR 00~2 022 

DEFSIOflflC DEFMRCRO 082 836 

DELF DEFSIOflflC 882 856 

DESIG-HRSH EXPR 803 826 

EFFECT DEFMRCRO 002 040 

GENRME EXPR 881 008 

HRSH EXPR 083 005 



HASH-RRRflY DEFMRCRO 003 002 

INCREMENT DEFSIDHRC 082 058 

INTEGERS-BETUEEN OEFTRIL 001 826 

LEXICAL -COMP ... EXPR .... 885 832 

LISTP DEFMRCRO 801 813 

LOGXOR DEFMRCRO 801 816 

MAKE-WRRP-STRCK EXPR 004 821 

MERGE DEFTfllL 881 819 

MUS2 EXPR .... 804 828 

NAMCOPY EXPR .... 801 002 

NAME-TYPE EXPR 001 885 

NCDRS DEFTfllL 805 829 

NODE -HASH EXPR 803 819 



NUM-COMP DEFMRC .. 

POP-QUERY EXPR 

PUSH-QUERY EXPR 

RANHASH DEFTfllL 

SCOHP EXPR 

SETF DEFMRCRO 

SUPER IOR-CLRSS OEFMRCRO 

SYMBOL -HRSH .... EXPR 

TRBLE-RSSCR EXPR 

VIEU-QUERY EXPR .... 

U-POP DEFSIOflflC 

U-PUSH DEFSIOflflC 

U-TOP DEFflflCRO 
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